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CHAPTER  I 


INTRODUCTION 

The  purpose  of  this  research  is  to  examine  the  theoretical 
possibility  of  using  pulsed  lasers  to  determine  the  velocity  struc- 
ture of  the  turbulence  associated  with  tornadoes^  particularly, 
the  funnel  cloud.  ^ p C 

The  basic  structure  of  a tornado  consists  of  a cloud  base, 
a funnel,  and  a debris  cloud.  The  cloud  base  may  be  visualized  as 
the  top  of  the  tornado,  connecting  the  clouds  to  the  tornado  funnel. 
The  funnel  extends  from  this  cloud  base  towards  the  ground;  if  it 
does  not  contact  the  ground  it  is  usually  referred  to  as  a funnel 
cloud.  If  the  funnel  does  come  in  contact  with  the  ground,  the 
severe  nature  of  its  turbulence  creates  a debris  cloud  which  may  ex- 
tend back  into  the  air  for  several  hundred  feet.  Physically,  the 
tornado  may  extend  from  a narrow  funnel,  with  a diameter  of  several 
meters,  to  much  larger  funnels  whose  damage  paths  indicate  diameters 
of  over  a kilometer.  Internal  velocities  of  particles  and  debris 
are  of  the  order  of  50-80  m/sec,  but  may  exceed  100  m/sec  (220  miles/ 
hour) . ^ Estimates  of  tornado  velocities  and  similar  turbulence  have 

^Robert  Davies-Jones  and  Edwin  Kessler.  "Tornadoes,"  Weather 
and  Climate  Modificxztion,  ed.,  Wilmot  N.  Hess  (New  York:  John  Wiley 
and  Sons,  Inc.,  1974),  p.  553. 
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come  from  direct  measurements  supplied  by  one  of  the  following 
methods:  doppler  radar,  doppler  lidart  aircraft  probe,  and  photo- 
gramraetry.  Each  of  these  techniques  is  limited  by  different 
factors:  The  radar  is  limited  by  a range  resolution  of  about  150  m. 
The  cw-lidar  used  by  Schwiesow,  et  at,,  is  limited  in  range  to 
about  a kilometer  due  to  power  and  mobility  requirements.^  Air- 
craft probe  is  considered  too  dangerous  for  tornadoes  and  photograra- 
metry  is  limited  dimensionally. 

Of  the  four  techniques,  the  lidar  appears  the  most  pro- 
mising. However,  the  range  limitation  remains  a serious  problem 
which  can  be  alleviated  by  using  a pulsed  laser;  a pulsed  laser 
has  the  capability  to  provide  considerably  more  pov.'er  than  a con- 
tinuous wave  laser  (cw-lidar]  and  should  thus  provide  more  range. 

Proposed  Scheme 

The  proposed  scheme  consists  of  analyzing  the  extinction 
and  backscatter  properties  of  mathematical  models  of  a tornado 
(based  on  probable  composition)  at  various  wavelengths  ranging  from 
0.5  ym  to  10.6  yra.  The  purposes  of  the  analysis  are  1)  to  deter- 
mine the  power  required  to  project  a beam  of  pulsed  laser  light 
through  the  intervening  atmosphere  to  the  tornado,  scatter  from  the 

^Robert  F.  Abbey,  Jr.,  Nineth  Conference  cm  Severe  Local  Storms 
(Norman,  Oklahoma,  1975),  p.  371-372. 

iHd.j  371. 
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particles  within  the  tornado,  and  return  to  a coaxially  mounted 
detector  with  enough  detectable  signal  to  allow  synthesis  of  the 
frequency  spectra,  2)  to  determine  which  wavelength (s)  is (are) 
the  most  conducive  to  this  type  of  analysis,  and  3)  to  determine 
the  optimum  pulse  width  required  to  obtain  the  desired  spectra. 

The  resultant  spectra  may  then  be  analyzed  to  obtain  the  Doppler 
shifts  and  thence  the  velocities  of  the  turbulence. 

Basic  Assumptions 

The  following  assumptions  served  as  a foundation  for  the 
calculations  and  results  of  this  investigation. 

1.  The  tornado  emanates  from  a cloud  structure  and,  until 
it  contacts  the  ground,  it  will  maintain  the  particulate  composition 
of  the  generating  cloud. 

2.  The  cloud  composition  consists  of  particles  ranging  in 
size  from  less  than  1 pm  to  about  25  urn  with  mode  radius  of  the 
order  of  several  micrometers  and  concentrations  which  depend  upon 
the  type  of  cloud  considered. 

3.  The  three  models  used  as  test  distributions  approximate 
extreme  conditions  of  cloud  composition  (relative  to  tornado-t>'pe 
turbulence) . 

4.  The  particles  within  the  cloud  are  spherical  and  par- 
tially absorbing.  (And  thus  the  scattering  which  predominates  is 
commonly  referred  to  as  Mie  scattering.) 

5.  The  scattering  that  occurs  from  a particular  particle 
is  independent  of  the  scattering  from  other  particles  in  the  cloud. 


CHAPTI-R  II 
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THEORETICAL  BACKCROUNP 


Laser  Propagation  in  the  Atmosphere 

The  propagation  of  light  emitted  into  the  atmosphere  by  a 
laser  is  unique  in  that  it  is  highly  directional  and  very  nearly  mono 
chromatic.  In  addition,  a laser  beam  has  a high  degree  of  both  spat- 
ial and  longitudinal  coherence  (coherence  times  can  be  on  the  order 
of  10"^  seconds)  and  radiates  as  a gaussian-spheri ca 1 light  beam 
(spherical  waves  having  a gaussian  variation  in  amplitude  across  the 
wavefront).'*^  This  light  beam  may  be  unpolarized,  pl.ine-polari zed , 
or  el  liptical ly-jiolari  zed  as  desired. 

An  optical  wave,  propagating  in  the  z-direction  with  fre- 
quency, u),  may  be  represented  mathematically  as 

<>(x,y,z,t)  “ E' (x,y)e^  (volts/meter) 


where  E'(x,y)  is  the  complex  wave  amjilitude  whose  (normalized  for  unit 

power  flow)  form  is 

/'j  y 7'.  2.  2 7 

^)  ^ exp(-i  j } exp  {-  — (volts/meter) 


*A.  E.  Siegman,  An  Intt\)(lut'tion  to  Liaero  and  Mjarnn  (New  York: 
McGraw-Hill  Co.,  1971)  .100. 

2r.  S.  Longhurst,  atui  Thyniaal  (.^rd  ed., 

London:  Longman  Group  Limited,  U)?.'^)  pp.  111-11.’^. 
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with  R representing  the  radius  of  curvature  of  the  spherical  wave 

and  W the  beam  radius.  ^ I 

i 

I 

One  may  apply  far-field  diffraction  theory  to  this  wave 
and  obtain  expressions  for  the  far-field  spot  size,  radius  of  curva- 
ture and  the  diffraction  angle.  These  expressions  are,  respectively, 

wcz)  . j 

2 

R(Z)  = 2[1+(tiW(o)“/XZ)‘]  = Z (meters);  Z (l-2b) 

2 

0 = W(Z)/R(Z)  = A/itW(o)  (radians);  2»  (l-2c) 

where  these  expressions  have  been  derived  on  the  basis  of  a '1/e' 
criteria  and  will  thus  contain  about  8(»  per  cent  of  the  full  gaus- 
sian  beam.  Note  that  the  laser  beam  docs  not  lose  its  sphcrical- 
gaussian  nature  when  transmitted  through  an  optical  lens.'*’^ 

In  order  to  achieve  a significant  amount  of  backscattcr 
at  the  distances  involved,  one  must  have  an  appreciable  scattering 

area.  Being  cognizant  of  the  further  requirement  that  one  of  the  ; 

objects  of  this  study  is  to  localize  the  area  of  the  tornado  from 
which  we  wish  to  obtain  the  data,  we  arbitrarily  selected  a scat- 
tering  area  of  10  m*"  as  a starting  point.  Relating  this  to  the  beam 


^A.  E.  Siegman,  Art  Iniroducrf:  wn  to  (jasora  and  Maeera,  p.  ,^07. 
'*rbid.,  308-309,  315. 

^F.  P.  liagliano  and  1).  .1.  Zaleckas.  "Laser  Processing  Funda- 
mentals," Ltaeva  in  Induatvy.  Edited  by  S.vS.  Charschan.  New  York: 
van  Nostrand  Reinhold  Company,  1972. 
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radius,  we  have,  at  the  scattering  area,  a beam  radius  of  about 
1.78  m.  Using  the  above  expressions  we  obtain  a spot  size  at  the 
transmitter  on  the  order  of  2 cm  for  a beam  of  wavelength  10.6  pm. 
Using  a typical  beam  waist  (of  a 10.6  Pm  laser)  of  about  5.0  mm  we 
require  optical  magnification  of  about  4X  at  the  transmitter  to 
achieve  the  desired  scattering  area. 

The  intensity  of  our  optical  wave  is  given  by  the  expres- 
sion 

I = (>,^*  = E’(x,y)E'*(x,y)  = exp  {-  } (watts)  (1- 

vW 

However,  we  are  primarily  interested  in  the  ratio  of  the  emitted 
intensity  to  the  received  intensity.  This  expression  is  derived 
from  consideration  of  the  elements  involved;  those  dictated  by 
the  field  of  view  of  the  detector,  the  medium  through  which  the 
beam  travels,  and  the  scattering  volume  itself.  For  a diffraction- 
limited  detector,  the  field  of  view,  q,  is  given  by 

a = X^/A^  (steradians)  (1- 

where  A is  the  wavelength  of  the  beam  and  is  the  capture  area  of 
the  detector. 

The  medium  in  this  case  is  the  atmosphere  between  the 
transceiver  and  the  tornado  and  it  is  characterized  by  a transmis- 
sivity, which  is  dependent  upon  the  constituents  of  the  atmos- 
phere. This  transmissivity  will  be  discussed  in  detail  shortly, 
but  for  our  purpose  here  we  shall  consider  it  as  a multiplicative 
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constant.  The  scattering  volume  is  accounted  for  in  the  following  t 

manner:  The  back-scattered  radiation  is  dependent  on  the  physical 

dimensions  of  the  scattering  volume,  the  transmissivity,  t^,  and 

the  volume  back-scattering  function,  of  the  volume  irradiated. 

2 

We  have  selected  a scattering  area.  A,  of  10m  . The  length  of  the 
scattering  volume,  is  in  principle,  a function  of  the  width  of 
the  pulse.  The  transmissivity  is  again  dependent  upon  the  consti- 
tuents cf  the  scattering  volume  and  is  given  by  the  expression 

= exp{-3g^j.L}  (dimensionless)  (1-5) 

where  6 ^ is  the  extinction  coefficient  derived  from  the  Mie 

ext 

theory  and  L is  as  defined  above.  The  volume  back-scattering  func- 
tion, B(n),  is  also  computed  from  the  Mie  theory.  Details  for  de- 
riving B(tt)  and  3^^^.  are  given  later. 

We  can  summarize  the  above  discussion  with  the  following 
relation,  noting  that  the  radiation  must  traverse  the  path  in  both 
directions, 

^RECEIVER  " H(X)ALnT^^3(ir)  exp{-2LBg^^)  (watts)  (1-6) 

2 

where  H(A)  is  the  spectral  irradiance  in  watts/m  and  the  other 
quantities  are  as  discussed  above. ^ 

The  field  of  view  for  a diffraction-limited  detector  is 
given  in  Table  I for  appropriate  values  of  wavelength  using  a 15.24 
cm  diameter  lens  for  the  aperature. 

®D.  B.  Rensch  and  R.  K.  Long.  "Comparative  Studies  of  Extinc- 
tion and  Backscattering  by  Aerosols,  Fog,  and  Rain  at  10.6  and  0.63  u", 

/pplied  Qptias  9,  No.  7 (July,  1970)  p.  1563. 
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Table  1.  Field  of  View  for  a 15.24  cm.  Diffraction 
Limited  Detector  at  Various  Wavelengths. 


Wavelength 

(Micrometers) 

Field  of  View 
(Steradians) 

10.6 

0.616x10'® 

4.0 

0.877x10'^° 

1.06 

0.616x10'^^ 

0.5 

0.137x10'^^ 
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Atmospheric  Attenuation 

Given  this  pulsed  laser  signal  emanating  from  the  trans- 
mitter into  the  atmosphere  towards  the  tornado,  we  seek  to  derive 
an  expression  for  the  attenuation  of  the  signal  by  the  intervening 
medium  and  the  turbulence  itself.  This  attenuation  is  a function 
of  the  nature  of  the  medium;  the  molecular  constituents  and  aero- 
sols (dust,  smoke,  carbon  particles,  etc.)  present  under  any  given 
circumstance.  For  most  conditions,  the  attenuation  assumes  an  ex- 
ponential form  given  by  an  extension  of  the  Lambert-Beers  law: 

= expC-oL)  (dimensionless) 

where  t is  the  atmospheric  transmissivity,  a is  the  attenuation  or 
A 

extinction  coefficient,  and  L is  the  path  length.  In  the  type  of 
media  commonly  referred  to  as  normal  or  clear  air  atmosphere,  the 
attenuation  coefficient  is  composed  of  four  parts,  i.e., 

a = o +K  +a. +IC.  (length)  ^ 

where  o,  and  k denote  scattering  and  absorption,  respectively.  The 
subscripts  a and  b refer  to  aerosol  and  molecular  effects,  respec- 
tively.® Of  these  four  parts,  only  two  have  significant  bearing 

for  the  wavelengths  we  are  considering  in  this  research;  these  are 
molecular  absorption  and  aerosol  scattering.  Molecular  scattering 
and  aerosol  absorption  prove  to  be  negligible  in  comparison. 


®Paul  W.  Kruse,  Laurence  0.  McGlauchl in , and  Richard  R.  McQuistan. 
Elemente  of  Infrared  Technology  (New  York:  John  Wiley  & Sons,  Inc., 
1962),  p.  162-192. 
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For  molecular  absorption,  the  components  of  consequence 
are  generally  considered  to  be  water  vapor  (H2O) , carbon  dioxide 
(C02)»  and  ozone,  but  the  latter  is  dropped  for  two  reasons:  the 
ozone  concentration  at  the  altitudes  we  are  interested  in  (sea 
level  to  12  kilometers)  is  negligible  and  the  absorption  coeffic- 
ient for  wavelengths  of  0.90  pm  and  above  is  zero.^>®  N2  and  O2, 
which  together  account  for  about  98  per  cent  of  the  volume,  do  not 
exhibit  absorption  bands  since,  being  homonuclear,  they  do  not  pos- 
sess either  a permanent  or  induced  electric  moment.  To  minimize  the 
effect  of  molecular  absorption,  one  must  avoid  the  sharp  absorption 
bands  associated  with  water  vapor  (1.1  pm,  1.38  pm,  1.87  pm,  2.7  pm, 
and  6.0  pm)  and  with  CO2  (2.7  pm,  4.3  pm,  and  14.5  pm).  If  these 
bands  are  avoided,  the  transmissivity  will  be  significantly  greater.^ 
Molecular  scattering  (Rayleigh  scattering)  is  negligible  in 
our  case  as  may  be  seen  from  the  following  considerations:  the 
molecular  diameters  of  the  most  prevalent  constituents  of  the  atmo- 
sphere are  all  on  the  order  of  a few  angstroms  (10  meters). 

The  attenuation  coefficient  may  be  ascertained  by  the  following 
formula: 

^William  K.  Pratt.  Laser  Communication  Systems  (New  York:  John 
Wiley  § Sons,  Inc.,  1969),  p.  129. 

^Louis  Alterman  and  Robert  B.  Tollin.  "Atmospheric  Optics," 
Bankbook  of  Geophysics  and  Space  Environments,  Ed.  Shea  L.  Valley 
(New  York:  McGraw-Hill  Co.,  1965),  p.  7-30. 

^Kruse,  Elements  of  Infrared  Technology,  pp.  163-164. 

iORobert  C.  Weast.  Ed.,  Handbook  of  Chemistry  and  Physics  (54th 
Ed.,;  Cleveland:  CRC  Press,  1973),  p.  F191. 
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Ojj  = 0.827  N Ap^A-**  (length)"^  . (2-1) 

where  n is  the  number  of  particles  per  unit  volume  in  the  path, 

Ap  is  the  cross  sectional  area  of  the  scatter  (a  fraction  of  the 
circular  cross-section) , and  A is  the  wavelength  of  the  incident 
light. ^ This  formula  yields  an  attenuation  coefficient  on  the 
order  of  10"^  for  circular  cross  sections.  The  listed  value  of 
the  attenuation  coefficient  for  0.5  urn  is  1.716 x 10"^  km’^;  in 
comparison,  the  listed  value  of  attenuation  coefficient  for 
aerosols  at  0.5  um  is  0.167  km' ^ and  therefore  the  molecular 

g 

scattering  is  about  10  per  cent  at  this  wavelength.  As  wave- 
length increases,  decreases  as  the  fourth  power  while  de- 
creases only  by  an  order  of  magnitude  (as  will  be  shown  directly) 
in  the  range  from  0.50  urn  to  10.6  urn.  Thus  we  neglect  the  atten- 
uation due  to  molecular  scattering  at  these  wavelengths  in  the 
normal,  clear  air  atmosphere. 

Aerosol  attenuation  in  the  atmosphere  may  be  described 
by  the  Mie  theory.  Due  to  the  particular  nature  of  the  aerosols, 
aerosol  absorption  is  again  neglected  (as  it  is  small  in  compari- 
son 


^Pratt,  Laaer  Corrmmioation  Systems,  p.  131. 

®Elterman  and  Toll in.  Handbook  of  Geophysics  and  Space 
Envivorments , pp.  7-30. 
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tothe  effect  of  scattering)  and  an  empirical  relationship,  based 
on  this  theory,  has  been  derived  which  gives  good  results  for 
clear  air  scattering.  This  relationship  is: 


a 


a 


5.91  , X .-O.SSSV^^^ 
V ^0.55^ 


(kilometers)  ^ (2-2) 


where  V is  the  visual  range  in  kilometers ^ However,  this  formula 
should  be  regarded  as  a minimum  since  it  does  not  include  variations 
in  local  meteorological  conditions.  A review  of  representative  values 
for  the  aerosol  attenuation  coefficient  indicates  variation  in  values 
of  almost  twice  that  given  by  Equation  (2-2)  for  1.06  pm,  over  five 

1 2 

times  greater  for  4 pm  and  almost  thirteen  times  greater  for  10.6  pm. 
Table  2 reflects  these  variations  for  a visibility  of  10  km. 
Similarly,  local  meteorological  conditions  affect  molecular  absorp- 
tion. A variation  in  the  atmospheric  absorption  of  CO2  of  0.08db/km 

between  January  and  July  has  been  reported  by  Yin  and  Long  for  a 

1 3 

wavelength  of  10.6  pm.  Water  vapor  absorption  is  a function  of 
the  water-vapor  concentration  which  in  turn  depends  on  both  relative 
humidity  and  temperature.  The  figures  for  molecular  absorption  in 
Table  2 have  been  obtained  using  concentrations  of  CO2  and  water 
vapor  of  700  atm-cm  and  15.0  precipitable-cm  respectively,  in 

Kruse,  Elements  of  Infrared  Technology,  p.  191. 

12 

Douglas  P.  Woodman.  "Limitations  in  Using  Atmospheric  Models 
for  Laser  Transmission  Estimates."  Applied  Optics  12:10,  2193 
(October,  1974). 

1 3 

P.K.L.  Yin  and  R.  K.  Long.  "Atmospheric  Absorption  at  the 
Line  Center  of  P(20)  CO-  Laser  Radiation."  Applied  Optics  7:8,  1551 
(August,  1968). 
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the  transmission  tables  of  Wyatt,  Stull  and  Plass . The  CO2 

concentration,  atmosphere-centimeter  is  defined  as  "...  the  length 
in  centimeters  of  a column  whose  volume  would  contain  the  same  number 
of  molecules  of  CO2  gas  at  NTP  as  a column  of  the  same  cross  section, 
but  of  arbitrary  length,  pressure,  and  temperature."  The  H2O  concen- 
tration, precipitable-centimeters , is  defined  as  "...the  thickness 
of  the  liquid  that  would  be  formed  if  all  the  water  vapor  traversed 
by  a beam  of  uniform  cross  section  were  completely  condensed  in  a 
container  of  cross  section  equal  to  that  of  the  beam."'**  The  700  atm 
-cm/km  gives  a CO2  absorption  coefficient  equivalent  to  0.3  db/km, 
which  is  consistent  with  the  data  of  Yin  and  Long.  The  15.0  preci- 
pitable- cm/km  corresponds  to  the  data  achieved  by  McCoy,  Rensch, 
and  Long.^^  The  values  in  Table  2 seem  to  cover  most  of  the  published 
values  for  the  effects  noted;  however,  it  is  stressed  that  local 
meteorological  conditions  will  change  these  figures,  especially  those 
of  water  vapor  absorption. 


l'*John  N.  Howard,  John  S.  Caring,  and  Russell  G.  Walker. 
"Transmission  and  Detection  of  Infrared  Radiation."  Handbook  of 
Geophysics  and  Space  Environments.  Edited  by  Shea  L.  Valley,  New 
York:  McGraw-Hill  Co.,  1965. 

^^Philip  J.  Wyatt,  Robert  V.  Stull, and  Gilbert  N.  Plass. 

"The  Infrared  Transmittance  of  Water  Vapor."  Applied  Optics  2:2, 

229  (February,  1964) . 

^^Philip  J.  Wyatt,  Robert  V.  Stull,  and  Gilbert  Plass. 

"The  Infrared  Tramsmittance  of  Carbon  Dioxide."  Applied  Optics  3:2, 
243  (February,  1964). 

^^John  H.  McCoy,  David  B.  Rensch,. and  Ronald  K.  Long.  "Water 
Vapor  Continuum  Absorption  of  Carbon  Dioxide  Laser  Radiation  Near 
lOp."  Applied  Optics  8:7,  1471  (July,  1969). 
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Table  2.  Transmission  Loss  of  the  Intervening  Atmosphere. 


Wavelength 

(micro- 

meters) 

Aerosol 

Scattering 

(db/km) 

Molecular 

Absorption 

(db/km) 

CO2  H^O 

Transmission 

Loss 

(db/km) 

10.6 

1 

o 

0.3 

1.6 

1.94-2.39 

4.0 

.14-. 75 

2.16 

1.0 

3.30-3.91 

1.06 

.74-1.44 

0.0 

0.09 

0.83-1.53 

0.5 

1.91-2.14 

0.0 

0.0 

1.91-2.14 
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Mie  Scattering  Theory 

Scattering  theory  has  developed  as  a natural  outgrowth  of 
the  study  of  light.  As  such,  it  was  greatly  influenced  by  Maxwell's 
electromagnetic  theory  of  light.  Mie  solved  Maxwell's  equations 
for  various  types  of  homogeneous  spheres  using  approximate  boundary 
conditions.  The  full  details  may  be  found  in  the  original  paper 
published  in  1908  {Ann.  Physik  26j  377).  A concise  outline  is  found 
in  many  books  including  van  de  Hulst.  The  Mie  solution  is  considered 
the  analytical  solution  as  it  is  a formal,  mathematical  analysis  of 
scattering  from  a sphere.  It  is  applicable  to  all  sizes  of  spheres 
and  uses  a size  parameter,  x = kr  = 2iir/X,  where  r is  the  radius  of 
the  sphere  and  X the  wavelength  of  the  incident  light.  Several  ap- 
proximations have  been  developed  including  the  two  extremum  approxi- 
mations, the  Rayleigh  (scattering)  approximation  (x<<l)  and  geome- 
trical optics  (x^®) . The  Mie  solutions  confirm  these  approximations 
which  in  turn,  are  normally  used  in  lieu  of  the  Mie  solutions  simply 
because  of  the  computational  analysis  involved  with  the  latter.  How- 
ever, no  acceptable  approximation  has  been  developed  for  the  size 
range  we  are  investigating  (x=l)  so  we  adopt  the  Mie  solution  and  proceed. 

The  results  of  a single,  independent  scattering  process  of 
a scalar  wave  may  be  generally  expressed  in  the  following  manner: 


[^2  ^31 

-ikr+ikz 

fE„1 

y = 

6 

ikr 

El. 

(volts/meter) 

o 

(2-3) 

where  the  SL  are  amplitude  functions  dependent  upon  9 and  <f>,  the  angle 
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between  the  directions  of  propagation  of  the  incident  light  and 
the  scattered  light  and  the  azimuthal  angle  respectively;  E,,  and 
are  the  electric  field  components  parallel  to  and  perpendicular  to 
the  direction  of  propagation,  z;  t is  the  radial  vector  perpendicu- 
lar to  z;  and  k is  the  propagation  constant.  For  spherical  scat- 
terers,  and  are  zero,  so  one  of  our  primary  objectives  is  to 
obtain  and  S^. 

They  have  been  derived  by  applying  boundary  conditions  to 
the  general  solutions  of  the  scalar  wave  equation  in  polar  coordi- 
nates. Noting  that  our  interest  in  the  scattered  wave  is  in  the 
far  field,  approximations  for  the  spherical  Bessel  functions  are 

applicable.  The  resulting  normalized  components  of  the  electric 
and  magnetic  fields,  sve  (using  polar  notation) 

E„  = H . = - e (coS()))S.,(e)  (dimensionless)  (2-4a) 

o ^ Kr  z 

-E^  ~ ^0  ~ ® ^*^^^^***^(sin()))Sj^(0)  (dimensionless)  (2-4b) 

oo  2ll^  1 

^1^®^"^  nfn^  (dimensionless)  (2-5) 

n=l  '■  ^ 

oo  2 1 

^ n('n+l)  (dimensionless)  (2-6) 

n-1  ^ 

where  u = circular  frequency  (radians/sec), 

“ ih  Pit™*®)  . 

TjjCcose)  ■ ^ p^(cose)  , 
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dependent  upon  the  complex  index  of  refraction,  m,  and  the  size 
parameter,  x.^®  They  may  be  expressed  most  compactly  in  terms  of 
the  Ricatti-Bessel  functions,  are, 

Ajj(y)'J'j^(x)-mr  (X) 

a * (dimensionless) 

" A^(y)Cj^(x)-mC;j^(x) 


b = 
n 


mA  (y)’f  (x)-'!"  (x) 
n n n 

«AnCy)c„(*)-i;;(x) 


(dimensionless) 


where. 


y = mx 


■ y (y) 


n 


’i'"’  ■ ii  ’n'*' 


K'  W = ^ (x) 


= Jf 

5„fx)  -/f  (x)l 


(all  dimensionless) 


C.  van  de  Hulst.  Light  Scattering  by  Small  Particles. 
New  York:  John  Wiley  5 Sons,  Inc.,  1957. 

Deirraendj ian.  Electromagnetic  Scattering  on  Spherical 
Poly dispersions.  New  York:  American  Elsevier  Publishing  Company, 
Inc.,  1969. 
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At  this  point,  we  can  determine  the  nature  of  the  electro- 
magnetic fields  of  an  incident  beam  scattered  from  a homogeneous 
sphere  of  radius,  r,  at  an  angle,  0,  from  the  incident  beam  and  at 
an  azimuthal  angle,  Furthermore,  with  the  additional  assumptions 
of  single  and  independent  scattering,  we  can  apply  the  results  to  N 
particles  of  radius,  r,  by  simply  invoking  the  principle  of  super- 
position and  adding  the  contributions  from  each  particle.  Maintaining 
these  assumptions,  we  can  also  extend  the  theory  to  include  a distri- 
bution of  particles  by  computing  the  Mie  amplitude  factors  and  coef- 
ficients of  each  representative  particle  and  summing  the  resultant 
contributions  weighted  by  their  relative  occurrence  in  the  distribu- 
tion. In  these  cases,  we  must  extend  the  notation  of  the  S^(0)  to 
include  the  functions  of  x and  m,  S^(0)  becomes  S^(x,m,0). 

There  is  a further  aspect  we  have  not  addressed  explicitly 
as  yet  and  that  is,  the  extent  to  which  the  incident  light  is  scat- 
tered and/or  absorbed.  As  intimated  earlier  in  the  case  of  molecu- 
lar scattering,  only  a fraction  of  the  incident  light  is  scattered 
and/or  absorbed.  The  fraction  that  is  scattered  (absorbed)  is 
equivalent  to  the  scattering  (absorption)  cross  section  divided  by 
the  geometrical  cross  section.  It  is  referred  to  as  the  scattering 
(absorption)  efficiency  factor,  » (the  K is  replaced  by  a 

Q in  van  de  Hulst's  notation).  Another  convenient  parameter  is  the 
total  extinction  efficiency  factor,  which  is  the  algebraic  sum 

of  the  scattering  and  absorption  efficiencies.  and  niay  he 

derived  from  the  integration  of  the  Mie  amplitude  functions  Equations 
(2-5)  and  (2-6)  and  are  given  by 
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I 


K = -^  I (2n+l)(|a  |2+|b  |2)  (dimensionless)  (2-7) 

SC  1 n n 

* n=l 

00 

K ^(x,m)  = ^ 7 (2n+l)Re{a  +b  } . (dimensionless)  (2-8) 

ext  1 n n 

X H"  1 

Equations  (2-5)  through  (2-8)  are  considered  the  fundamen- 
tal parameters  of  the  Mie  Theory,  as  all  the  other  parameters  re- 
quired to  describe  the  intensity  and  polarization  of  scattered 
light  may  be  derived  from  these  four  parameters .2° 

We  are  now  in  a position  to  describe  those  functions  pre- 
viously referred  to  as  the  extinction  coefficient,  ft  . and  the 

ext 

volume  backscattering  function,  Biir)-  These  functions  are  deter- 
mined from  the  particle  size  distribution,  n(r),  the  scattering  ef- 
ficiency and  as  yet  undefined  quantities  known  as  the  Mie  intensity 
functions,  ij , which  are  derived  from  the  amplitude  functions  Equa- 
tions (2-5)  and  (2-6)  by  means  of  the  traditional  method  of  obtaining 
intensities  from  amplitudes,  that  is 

ij(x,ra,0)  = SjS^*  (dimensionless) 

i2(x,m,0)  = (dimensionless) 

For  a single  particle,  is  just  the  extinction  cross  section, 

but  for  a distribution  of  particles,  0^^^  must  include  the  contri- 
butions from  all  particles.  Therefore  we  once  again  apply  the 
principle  of  superposition  and,  expressing  0^^^^  in  terms  of  the 
parameter  x, 


p. 
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X 

t max 

®ext  ” nCx)x2K^^^(x,m)dx  (length) 


-1 


(2-9) 


'X  . 

min 


where  C is  the  constant  obtained  when  converting  n(r)  and  r to  n(x) 
and  X respectively. 

Similarly,  B(ir)  is  found  to  be 
, X 

^ f _1  -1 

= I d|  n(x)i . (x,in,Tr)dx  (length)  ^(steradian) 

j=l  Jx  . ^ 

where  D is  the  constant  which  includes  the  conversion  factors.  This 
can  be  further  simplified  by  noting  that  i^  = i2  for  0 * tt  so  that 


X 

r max 


e(v)  = D' 


X . 

roin 


n(x)  i^(x,ra,Tr)dx  (length)  ^(steradian)  ^.^^(2-10) 


We  now  proceed  to  describe  the  distribution  function  n(r) 
which  specifies  concentration,  mode  radius  and  shape  of  the  curve 
depicting  the  particles  in  our  scattering  volume. 

The  Distribution  Function 

Precipitating  clouds  have  characteristics  which  vary  widely 
depending  on  the  location,  humidity  and  other  conditions.  Our  pro- 
cedure was  to  select  clouds  associated  with  tornadoes  in  the  local 
area  (Oklahoma),  determine  the  characteristics  of  these  clouds,  and 

P.  14. 

pp.  89,  119. 
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then  proceed  to  obtain  distribution  functions  which  closely  match 
those  characteristics. 

The  local  continental  clouds  are  believed  to  have  the  fol- 
lowing characteristics:  (1)  The  concentration  in  the  funnel  cloud 
has  a maximum  near  1000  particles/cm 3 and  a minimum  of  about  150 
particles /cm  3.  (2)  The  mode  radii  associated  with  the  above  con- 

centrations are  5 pm  and  15  pm  respectively.  (3)  The  relative  dis- 
persion, = (the  standard  deviation  of  particle  radii)/ (mode 
particle  radius),  on  the  order  of  0.366.  And,  (4)  The  significant 
portion  of  the  particles  being  in  the  range  of  1 pm  to  25  pm  in 
radius.  2 

The  distribution  function  used  was  first  proposed  by 
Deirmendjian  and  has  since  been  used  by  several  investigators  in 
similar  applications. Selected  on  the  basis  of  its  versa- 
tility and  relation  to  physical  interpretation,  that  function  is, 

n(r)  = ar“exp(-br^)  cm"3/pm  0<r_<oo  (2-11) 


where  the  four  constants  a,  b,  a,  y are  real  and  positive,  and  a is 


23joanne  Simpson  and  Victor  Wiggert.  "Models  of  Precipi- 
tating Cumulus  Towers."  Monthly  Weather  Review  97:7,  471  (July, 
1969). 

John  McCarthy.  Personal  Interview.  September,  1975. 

D.  Deirmendjian.  Scattering  and  Polarisation  Properties  of 
Polydispersed  Suspensions  with  Partial  Absorption.  Report  prepared 
for  United  States  Air  Force  Project  Rand,  Memo#  RM-3228-PR,  July, 
1962.  Santa  Monica,  CA.'Rand  Corporation,  1962. 

^®R.  C.  Anderson  and  E.  W.  Browell.  "First-  and  Second-Order 
Backscattering  from  Clouds  Illuminated  by  Finite  Beam."  Applied 
Optics  11:6,  1345  (June,  1972). 

D.  B.  Rensch  and  R.  K.  Long.  "Comparative  Studies  of  Ex- 
tinction and  Backscattering  by  Aerosols,  Fog,  and  Rain  at  10. 6p  and 
0.63p."  Applied  Optics  9:7,  1563  (July,  1970). 
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an  integer.  Further  properties  are  noted  by  Oeirmendjian  in  the 
original  work  and  in  his  book  (Deirmendj ian,  1969);  however,  only 
the  values  of  the  constants  concern  us  here,  a and  y are  obtained, 
in  our  case,  by  trial  and  error  to  fit  the  relative  dispersion,  Dj,, 
and  the  estimated  range  of  the  particles,  "b"  is  obtained  by  dif- 
ferentiating n(r)  with  respect  to  r and  determining  the  absolute 
maximum;  it  is  given  by, 

b = (pm)'^  (2-12) 

where  r^.  is  the  mode  radius  of  the  distribution,  "a"  is  obtained 
by  integrating  n(r)  over  the  range  of  r and  solving  for  a particu- 
lar value  of  N,  the  concentration  of  particles.  Thus,  a is  given 
by 

a = ^Y^OL+'l/'y)  (particles/cm3)(pm)'^“^^^  (2-13) 

where  r(x)  is  the  gamma  function.^® 

Doppler  Effect 

To  determine  the  velocity  of  the  particles  from  which  the 
incident  light  beam  is  scattered  one  applies  the  principle  of  the 
Doppler  effect.  To  be  precise,  one  would  use  the  precepts  of  the 
Theory  of  Special  Relativity  to  obtain  the  frequency  shift  and 
thence  the  velocity.  For  incident  light  with  frequency  f,  the  ex- 
pression for  the  Doppler-shifted  frequency  f is 


2®Deirmendjian,  Scattering  and  Polax>ization,  pp.  2-3. 


d 


f.  . Hz 


where  0 = velocity  of  particle/speed  of  light  = v/c  and  the  top 
sign  is  used  for  motion  away  from  the  observer  while  the  bottom 
sign  is  used  for  motion  toward  the  observer. However,  we  are 
dealing  with  velocities  on  the  order  of  100  m/sec  so  we  use  the 
first-order  approximations  of 


f f 

f'  = -i—  = — i-T—  Hz 

1±S  r±v/c 


we  must  also  account  for  the  direction  of  the  velocity  of  the  parti- 
cle in  relation  to  the  receiver  since  it  is  the  velocity  component 
in  the  direction  of  the  receiver  that  accounts  for  the  change  in  the 
frequency.  Then,  one  obtains 


(2-14) 


1 ± — COS0 
c 


where  e is  the  angle  between  the  direction  of  the  particle  and  the 


receiver  as  in  the  sketch  below. 


RECEIVER 


PARTICLE 


29  A.  P.  French.  Special  Relativity.  New  York:  W.  W.  Norton 
and  Co. , Inc. , 1968. 

30  F.  P.  Gagliano  and  U.  J.  Zaleckas.  "Laser  Processing 
Fundamentals."  Lasers  in  Industry.  Edited  by  S.  S.  Charschan. 
New  York:  Van  Nostrand  Reinhold  Company,  1972. 
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To  determine  the  velocity  resolution  of  the  laser,  we 
must  take  into  consideration  the  spectral  width  , Sg,  of  the  laser 
since  a laser  does  not  emit  radiation  of  a single  frequency.  We 
proceed  as  follows:  A reasonable  velocity  resolution  of  Im/sec 
is  consistent  with  the  total  velocities  within  a tornado.  Using 
Equation  (2-14)  and  assuming  a velocity  component,  v^(=  vcos0), 
approaching  the  receiver  we  obtain  for  v^, 

f f'-f 

Vq  = (1-  p-)c  = -p-  c ^ -p-  c (meter/sec) 


To  obtain  a 1 m/sec  velocity  resolution  we  must  have  for  a wave- 

, _ ^ 2vf  2Clm/sec)2.83xl0^3/sec  , 

length  of  10.6  pm,  S„  = = — = 1.887x10^  Hz. 

* ’Be  3xl0®m/sec 

One  similarly  obtains  maximum  spectral  widths  of  5x10^  Hz, 

1.887x10^  Hz,  4x10^  Hz  for  the  wavelengths  of  4 pm,  1.06  pm,  and 
0.5  pm,  respectively.  According  to  De  Nicola,  these  spectral  widths 
are  well  within  the  capabilities  of  today's  lasers. 

The  bandwidth  of  the  receiver  is  also  affected  by  the  veloc- 
ity of  the  particle;  in  this  case  it  is  the  maximum  expected  veloc- 
ity of  lOOm/sec.  Applying  the  above  principles,  we  obtain  receiver 
bandwidths  of  1.887xl0^  Hz,  5x10^  Hz,  1.887xlo8  Hz,  and  4x10®  Hz  for 
wavelengths  of  10.6  pm,  4 pm,  1.06  pm,  and  0.5  pm,  respectively. 

A consequence  of  this  discussion  is  based  on  the  inherent 
property  of  the  Doppler  effect  to  account  for  the  velocity  component 
of  the  particle  parallel  to  the  direction  of  the  receiver.  Thus 
maximum  velocity  spectra  will  come  from  the  side  portions  of  the 
turbulence. 


0.  DeNicola.  "Reflection  and  Scattering"  Lasers  in  Indus- 
try. S.  S.  Charschan  (New  York:  Van  Nostrand  Reinhold  Co.,  1972). 
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Heterodyne  Detectors 

To  detect  the  backscattered  radiation,  we  propose  a coax- 
ially mounted,  heterodyne  detector.  A coaxial  mount  takes  full 
advantage  of  the  backscattered  radiation,  and  heterodyning  is  de- 
sirable for  a variety  of  reasons  which  include,  high  sensitivity, 
inherent  narrowband  detection  and  the  wide  range  and  availability 
of  solid  state  detectors  (with  their  higher  quantum  efficiencies) 
in  the  infrared  region. A further  advantage  of  heterodyning  is 
that,  as  will  be  shown,  one  can  effectively  eliminate  the  background, 
shot  and  thermal  noise  by  controlling  the  local  oscillator  power. 

The  essence  of  heterodyning  is  depicted  in  Figure  1.  Tlic 
signal  current,  i . , is  given  by  the  product  of  the  detector  pro- 

Slg 

portionality  constant,  3 = qn/hf^,  (the  responsivity)  the  gain,  G, 
and  the  total  signal  squared  (assuming  square  law  detection).  The 
result  is,  assuming  Aw<<a)^,  and  that  the  detector  is  unable  to 
follow  the  fast  variations  in  frequency, 

Es^  E^2 

i . ■v  eG(-^^  ♦ ♦ E E cosAut)  (Amperes) 

Slg  ■“  2 2 so 


eliminating  the  dc  terms,  the  power  delivered  to  the  load  resistor, 
Rl*  is 


S 


l2.  R, 
Slg  L 


G2b2e  2e  2r 

s o L 
2 


(watts) 


or,  in  terms  of  the  average  power, 

S - 2G262p  P R, 
SOL 


(watts)  .33  (2-15) 


32s.  S.  Charschan.  "Detection  and  Measurement,"  Lasers  in  Indus- 
try (New  York;  Van  Nostrand  Reinhold  Co.,  1972). 

33pratt,  Laser  Cemruniaation  Systems,  p.  187. 


Local  Oscillator 
Eg  Cos  tOgt 


Figure  1.  Basic  Heterodyne  Detector 
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The  total  noise  delivered  to  the  load  resistor  primarily 
arises  from  five  sources  which  manifest  themselves  as  shot  noise 
and  thermal  noise.  The  shot  noise  is  comprised  of  the  following: 
the  dark  current,  1^,  and  the  currents  generated  by  fluctuations 
in  the  background  radiation,  the  local  oscillator  and  the  carrier 
itself.  Included  in  this  type  of  noise  may  be  the  generation-recom- 
bination noise  for  photovoltaic  and  photoconductive  detectors.  The 
thermal  noise  is  the  familiar  Johnson  or  Nyquist  thermal  noise.  The 
total  expression  for  noise  in  terms  of  the  average  power  is 


Njp  = 2G2q[B(P5+PQ+Pg)  + Ip]BjpRj^  + 4kTBjp  (watts) 


where  Pg,  Pq,  Pg  represent  the  signal,  oscillator,  and  background 

noise  respectively,  k is  Boltzmann's  constant,  T is  the  temperature 

34 

in  degrees  Kelvin,  and  Bjp  is  the  bandwidth  of  the  IF  filter.  As 
stated,  making  the  power  of  the  local  oscillator  large  with  respect 
to  the  other  parameters  tends  to  make  those  other  parameters  negli- 
gible and 

Njp  = 2G2q6PoR^Bjp  . (watts)  (2-16) 


Combining  Equations  (2-15)  and  (2-16)  we  obtain  for  the  signal  to  noise 
ratio  (SNR) 


2G2g2PsPoRL  g Pg  n P3X 

2G2qBP^R^Bjp 


(dimensionless) 


which  becomes,  for  a sine  wave  plus  narrow-band  Gaussian  noise  in- 
put, using  the  square  law  detector. 


187. 
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2hcB 


IF 


(dimensionless) 


(2-17) 


Now,  B^p>2B_  and  so,  to  accommodate  the  Doppler  shift  as  previously 
discussed,  the  SNR  will  be  slightly  lower.  If  horaodyne  detection 
is  used,  the  bandwidth  of  the  local  oscillator  will  have  to  be  ex- 
tended similarly  for  the  doppler  effect. 

As  a minimum,  the  SNR  should  equal  unity,  so  we  may  calcu- 
late the  required  signal  power  from  Equation  (2-17)  using  the  re- 
ceiver bandwidths  calculated  in  the  previous  section.  We  obtain, 
for  a quantum  efficiency,  n,  of  1,  and  a wavelength  of  10.6  ym, 

P^  = 2Bjphf  = 4B^hf  = 4(1.887x10^  sec’ 1)  (6.626x10' 3**  joule- 
sec)  (2. 83x10' 3 sec’')  = 1.41x10’'^  watts.  Similarly,  we  obtain 
minimum  signal  powers  of  9.94x10’'^  watts,  1.41x10’'°  watts,  and 
6.36x10*'°  watts  for  the  wavelengths  of  4.0  ym,  1.06  ym,  and  0.5 
ym,  respectively. 

Ross  notes  that  "...communications  would  be  difficult  at 

SNR  = 1 and  hence  at  least  a factor  of  ten  is  required  for  adequate 
36 

communication."  Therefore,  our  signal  powers  should  be  an  order 
of  magnitude  larger  than  the  minimum  computed  above. 


3 5j2)td.,  pp.  187-190. 

36  • 

Monte  Ross.  Laser  Receivers  (New  York:  John  Wiley  and  Sons, 
Inc.,  1966),  p.  327. 
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CHAPTER  III 


PROCEDURES  AND  CALCULATIONS 


Computer  Implementation 

The  use  of  a computer  to  develop  the  Mie  parameters  and 
analyze  the  scattering  of  light  from  a distribution  of  particles 
is  a necessity  due  to  the  infinite  series  involved.  There  are 
tables  available  that  list  the  Mie  parameters  [see  (Deirmendj ian, 

1963)  for  example  and  (van  de  Hulst,1957)  for  a representative  listing 
of  various  sources] . However,  the  tables  are  unsuitable  for  some 
applications  due  to  truncation  of  the  values  or  to  interval  spacing 
which  is  too  large  for  present  day  applications  as  noted  for  example, 
in  works  by  Dave.^»2 

The  Mie  parameters  are  quite  adaptable  to  computerization  as 
has  been  shown  by  Dave.^  We  used  the  adaptations  of  Deirmendj ian 


'j.  V.  Dave.  "Effect  of  Coarseness  of  the  Integration  Increment 
on  the  Calculation  of  the  Radiation  Scattered  by  Polydispersed  Aero- 
sols." Applied  Optics  8:6,  1161  (June,  1969). 

^J.  V.  Dave.  "Effect  of  Varying  Integration  Increment  on  the 
Computed  Polarization  Characteristics  of  the  Radiation  Scattered  by 
Polydispersed  Aerosols."  Applied  Optics  8:10,  2153  (October,  1969). 

3j.  V.  Dave.  "Subroutines  for  Computing  the  Parameters  of  the 
Electromagnetic  Radiation  Scattered  by  a Sphere."  Report  prepared 
for  IBM,  Rep.  320-3237.  Palo  Alto,  CA:  IBM  Scientific  Center,  1968. 
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(see  Appendix  A)  to  determine  the  fundamental  Mie  parameters  as 
specified  in  Equations  (2-5)  through  (2-8).  The  interval  spacing 
was  approximately  0.25  wra.  Using  this  interval,  we  proceeded  to 
compute  the  distribution  function,  n(x)  and  ultimately  the  extinc- 
tion coefficient,  6^^^.  and  the  volume  backscattering  function, 

6(n),  using  numerical  integration  (trapezoidal  rule)  for  all  the 
distributions.  The  program  then  computes  the  received  intensity 
(neglecting  atmospheric  transmission  effects)  as  given  by  Equation 
(1-6),  for  one-meter  increments  of  the  length  of  the  scattering  volume. 

Appendix  B has  a flowchart  of  the  computer  program.  It  is 
not  all  inclusive,  but  is  given  as  an  aid  in  deciphering  the  program, 
^pendix  C has  the  complete  computer  program  printout.  For  simplic- 
ity, it  is  written  in  the  Fortran  Language,  using  a WATFIV  compiler. 


Prograumning  the  Mie  Parameters 
The  calculations  proceed  in  the  following  manner:  After 
initialization  of  the  various  parameters  an  "x-Do  loop"  is  established 
for  each  required  x.  The  end  point  values  of  x and  the  step  size.  Ax, 
are  listed  in  Table  3 with  their  corresponding  values  of  r and  Ar  for 
each  wavelength  considered.  For  each  x,  the  initial  values  of  A^, 

are  computed  as  A^,  W^,  respectively  to  avoid  using  nega- 
tive and  zero  indices.  A nested  " -Do  loop"  is  then  created  to  com- 
pute the  Mie  functions  W , n , t , A , a and  b and  finally  the  Mie 

nnnnnn  ' 

parameters  Sj(it)  [=82(11)],  K^^(x,m),  values  of  n from 

3 (corresponds  to  n = 1 as  explained  above)  to  a value  n = 1.2  +12.** 

**D.  Deirmendjian  and  R.  J.  Clasen.  Light  Saatteving  on  Partially 
Absorbing  Homogeneous  Spheres  of  Finite  Size.  Report  prepared  for  United 
States  Air  Force  Project  Rand,  #R-393-PR,  February,  1962.  Santa  Monica, 
CA:  Rand  Corporation,  1962. 
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This  latter  value  is  an  Iteration  larger  than  that  recommended  by 
Deirmendjian  to  insure  proper  convergence  of  the  finite  series  (after 
raising  the  indices).  Larger  numbers  were  tried  but  failed  to  alter 
the  results,  thus  verifying  the  convergence.  Hie  "n-l)o  loop"  ends 
after  those  iterations  are  complete  for  the  named  parameters. 


Determining  the  Scattering  .and  nxtiiicti on  Properties 
After  computation  of  the  Mie  parameters,  the  extinction-, 
scattering-,  backscatter-coefficicnts  and  tlie  distribution  functions, 
were  computed  for  each  x in  intervals  of  Ax  according  to  tlie  following 
scheme:  The  form  used  for  the  extinction  coefficient  is 


r 

> max 

6 =•  n r^n(r)K  (m,x)dr  (length) 

ext  Jr 

min 


To  express  this  witii  one  variable  x,  the  variables  r'  , n(r),  dr  must 
be  converted  to  expressions  of  the  variable  x. 

r^  •=  (r)'^  * (micrometers)- 

K 4.  IT 

, , ,xX,  , . ,xX,Y, 

n(r)  » n(:T-)  • a(:;r  ) exp  -b(^)  ) 

4.1T  4.1T  4:11 

Jr  - — - dx  (micrometers) 


Therefore, 


r2n(r)dr  » (-^-) 


3*a 


ax 


dx 


and. 


( 


■■ 


m 
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W " j ax“^^exp[-b(^)\^]Kg^^(jn.x)dx  fmeter)"^ 

^in 

which  is  in  the  form  previously  noted  in  Equation  (2-9),  the  addi- 
tional constant  of  10“ ® occurs  due  to  changing  the  units  so  that 

is  expressed  in  (meters)"^.  Note  that  A is  expressed  in  microns. 
Similarly  we  find  the  expression  for  the  scattering  coefficient, 

B , to  be 
^sc 

*max 

Ssc  = ^ x“^^exp[-b(^)V]K^^(m,x)dx 

(meter) 

For  the  backscatter  coefficient,  we  proceed  along  the  same 
lines  by  starting  with 

S(it)  = IT  r2n(r)K^(m,x)dr  (length)  -(steradian) 

Ir  . 
min 

where  Kj^,  the  backscattering  efficiency,  is  defined  as 

A ^ 

= — )Sj(ir)|  (dimensionless).^ 


Then  we 

8(11)  = 


find,  converting  to  the  variable  x as  before,  that 

X 

. A 1 3-*-a  , , n- 6 A f Cl  ru/-^\^Yllcr-v|2j 

^2ii^  (10  *>)  ax  exp[-b(jjj-)  x']|S^(it)|  dx 

* X 

min 


(meter) 


(steradian) 


-1 


^D.  Dcirmendjian.  Electromagnetic  Scattering  on  Spherical 
Potydiepersione.  New  York:  American  Elsevier  Publishing  Company, 
Inc.,  1969. 
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The  "x  Do- loop"  continues  with  the  computation  of  the  sum- 
mation terms  for  integration  by  the  trapezoidal  rule,  which  may  be 
expressed  in  the  following  form: 

.b 

f(x)dx  = ~ [%f(a)  + f(xp+f(x2)  + ”--^f(x^_j)+%f(b)]  .6 

^ a 

The  program,  as  currently  established,  computes  the  integra- 
tion for  the  extinction  coefficient  (BEXT  1,2, 3, 4),  the  scattering 
coefficient  (BSC  1,2, 3, 4),  and  the  backscatter  coefficient  (BPI  1,2, 
3,4)  for  four  different  distribution  functions  and  two  different 
limits  (each  initiated  at  the  same  point,  but  having  different 
maximums,  LAST  1,2):  The  first  set  of  limits  was  used  for  verifi- 
cation against  known  results  for  one  of  the  distribution  functions 
and  the  second  set  of  limits  was  used  for  the  set  of  distribution 
functions  that  represent  the  probable  composition  of  the  tornado 
funnel.  In  its  present  state,  the  program  is  easily  rearranged  to 
compute  the  coefficients  mentioned  in  any  combination  of  zero  to 
four  distribution  functions  in  one  set  of  limits  with  the  remaining 
distribution  functions  in  the  second  set  of  limits.  In  addition, 
one  may  easily  expand  (or  contract)  either  the  number  of  distribu- 
tion functions  or  the  munbcr  of  sets  of  limits  or  the  number  of 
functions  obtained  by  integration  or  all  three.  This  is  done,  by 
adding  the  applicable  distribution  functions  (FUN  5,6,...)  and/or 
integration  loops  (BEXT  5,6,...;  BSC  5,6,...;  etc.),  and  then 

®E.  Kreyszig.  Advanced  Engineering  Mathematics.  3rd  Ed. 

New  York:  John  Wiley  and  Sons,  Inc.,  1972. 


adjusting  the  format  statements  for  readout,  to  meet  the  specific 
need.  The  "x-Do  loop"  is  completed  by  a printout  of  x and  its  as- 
sociated Mie  parameters  at  selected  intervals  of  x. 

Determining  the  Receiver  Intensity 
The  program  continues  by  computing  the  receiver  intensity 
(Backscatter) , as  given  by  Equation  (1-6)  except  for  the  transmis- 
sivity of  the  intervening  atmosphere,  for  1 -meter  increments  of  the 
length  of  the  scattering  volume  (Resolution)  for  each  distribution 
function,  given  a specific  field  of  view  and  the  area  of  the  scat- 
tering volume  (this  quantity  is  entered  as  the  single  entity  "OMEGA" 
in  the  program) . [Appropriate  controls  have  been  entered  here  to 
keep  from  exceeding  the  negative  exponential  limit  of  the  computer.] 
The  results  currently  programmed  for  printout  are: 

(1)  X and  its  associated  Mie  parameters  for  selected  intervals  and 

(2)  The  extinction-,  scattering-,  and  backscattering-coefficients  to- 
gether with  the  Resolution,  solid  angle  (field  of  view--scattering 
area  function,  "OMEGA")  and  intensity  ratio  (Backscatter),  for  each 
distribution  function  identified  with  an  index  of  J = 1,2,3,  or  4. 

The  Distribution  Functions 

An  integral  portion  of  the  calculations  are  the  distribution 
functions.  Four  functions  were  utilized;  the  first  (FUN  1)  for 
verification  of  the  program  against  results  obtained  by  Carrier,  et 
at.,  for  their  "STRATUS  I"  cloud  (the  closest  distribution  to  our 
other  functions  in  relative  shape,  concentration,  etc.).^ 

^L.  W.  Carrier,  G.  A.  Cato,  and  K.  J.  von  Essen.  "The  Back- 
scattering  and  Extinction  of  Visible  and  Infrared  Radiation  by  Sel- 
ected Major  Cloud  Models."  Applied  Optica  6:7,  1209  (July,  1967). 
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The  last  three  functions  represent  probable  distributions  of  parti- 
cles in  a tornado  funnel  (which  has  not  contacted  the  ground) . The 
functions  are  plotted  in  Figures  2 through  5 and  the  associated 
constants  for  the  distribution  function,  n(r)  = ar°exp[-brT]  are 
listed  in  Table  4. 

Figure  2 contains  two  plots,  one  of  the  original  distribu- 
tion used  by  CARRIER,  et  at.  and  the  other  the  distribution  func- 
tion we  used  to  approximate  the  original.  The  results  obtained, 
as  the  data  indicates,  are  very  close  to  the  values  obtained  ori- 
ginally for  all  values  of  wavelength  used  in  this  study.  Figures  5 
and  5 represent  the  probable  maximum  and  minimum  concentrations, 
respectively,  of  the  expected  distributions.  Figure  4 represents 
a median  concentration.  The  concentration,  mode  radii,  and  rela- 
tive dispersion  are  indicated  on  each  of  the  figures. 

Program  Verification 

During  the  initial  development  of  the  program,  the  funda- 
mental Mie  parameters  defined  earlier  and  described  by  Equations 
(2-5)  through  (2-8)  were  rigorously  checked  against  published  re- 
sults for  the  10  urn  range  using  an  index  of  refraction  of 
m = 1.212-i0.0601.®'^* The  values  computed  matched  the  values 
published  by  Deirmendjian  to  the  accuracy  quoted  and  agreed  with 
those  of  Gumprecht  and  Sliepcevich.  For  the  wavelengths  of  0.5  pm 

®D.  Deirmendjian.  Eleatvomagnetic  Saatteving,  p.  30. 

9D.  Deirmendjian.  Tables  of  Mie  Scattering  Cross  Sections  and 
Amplitudes.  Report  prepared  for  United  States  Air  Force  Project  Rand, 

#R-407-PR,  January,  1963.  Santa  Monica,  CA:  Rand  Corporation,  1963. 

^®R.  0.  Gumprecht  and  C.  M.  Sliepccvish.  Light-Scattering  Func- 
tions for  Spherical  Particles.  Willow  Run  Research  Center:  Univers- 
ity of  Michigan  Press,  1951. 
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Program 

Name 

Mode 

Radius 

(lira) 

a 

b 

a 

Y 

FUN  1 

3.5 

23.480232 

.42341 

3 

1 

FUN  2 

5 

1.0 

0.1 

5 

2 

FUN  3 

10 

7. 81248x10' 3 

0.025 

5 

2 

FUN  4 

15 

1.1098x10'^ 

y.goixio'** 

8 

3 

n(r)  Absolute  concentration,  number  -cm" 
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Figure  3.  First  probable  funnel  cloud,  distribution  #2 


^ ^ 
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n(r)  Absolute  concentration,  number  -cm'^-pm 
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and  1.06  vim  with  indices  of  refraction  of  1.34  and  1.33  respectively, 
the  computed  parameters  were  printed  out  at  intervals  to  preclude 
volumes  of  numbers.  In  all  cases,  the  values  printed  correspond 
to  the  published  results  of  Gumprecht  and  Sliepcevich  and  matched 
the  values  computed  from  the  exact  Mie  theory  as  cited  by  Houghton 
and  Chalker.^^  Published  results  for  index  of  refraction 
m = 1.353-i0.0059(A  =4.0  Mm)  were  not  available  for  comparison,  how- 
ever, the  computed  results  did  correspond  quite  well  in  degree  of 
magnitude  and  inflection  with  the  values  published  for  m = 1.33. 

The  values  obtained  for  the  extinction  and  backscatter  co- 
efficients were  checked  by  inserting  a test  distribution  which  ap- 
proximated the  "Stratus  I"  cloud  of  Carrier,  et  at.,  (See  Figure  2).^** 
The  values  obtained  are  listed  in  Table  5 along  with  the  values  ob- 
tained by  Carrier,  et  al.,  and  the  percentage  difference  computed 
using  the  latter  values  as  the  standard.  The  differences  in 
are  caused  by  the  narrower  relative  dispersion  at  the  median  of  the 
spectrum  of  Carrier,  which  tends  to  increase  3^^^. 

H.  G.  Houghton  and  W.  R.  Chalker.  "The  Scattering  Cross 
Section  of  Water  Drops  in  Air  for  Visible  Light,"  Joujcmal  of  the 
Optical  Society  of  America  29:11,  955  (November,  1949). 

^ ^Gnimprecht  and  Sliepcevich,  Light-Scattering  Functions,  p.  4. 

^^Houghton  and  Chalker,  "Scattering  Cross  Section"  pp.  955-957. 

^‘‘Carrier,  "Backscattering  and  Extinction"  p.  1212. 
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Table  5.  Comparison  of  Extinction  and  Scattering  Properties 
for  the  Test  Distribution. 


Wavelength  Computed 

(Micrometers)  Values 


Per 

Carrier,  et  dlJ  Cent 

Difference 


10.6 

®ext 

3.92x10~2 

m-  ^ 

3(w) 

6.41x10-5 

m"  ^ 

4.0 

^ext 

8.46x10-2 

m"  ^ 

3 (If) 

1.25x10-3 

ffl-  ^ 

1.06 

®ext 

6.80x10-2 

m-  ^ 

3 (if) 

3.59x10'3 

m"  ^ 

4.28x10-2 

m-^ 

-8.4 

sr-l 

7.42x10-5 

m-  ^ 

sr-  ^ 

-13.6 

9.01x10'2 

m’  ^ 

-6.1 

sr'^ 

1.47x10-3 

m- 1 

sr-l 

-15.0 

6.97x10-2 

m-  ^ 

-2.4 

sr~^ 

3.08x10-3 

m-^ 

sr- 1 

+ 16.6 

0.5 


®ext 

S(ir) 


6.69x10"2  m"l 
3.13xl0‘^m-^r-^ 


I 


6.54x10‘2  m’l 
3.49x10‘3  ra'^sr'l 


-2.2 

+ 11.5 


4S 


The  positive  increase  in  the  values  of  6 (if)  at  1.06  wn  and  0.5  urn 
is  predicted  by  the  analysis  of  Carrier,  and  is  presumably  due  to 
a combination  of  broader  spectrum  and  larger  particle  concentra- 
tion in  the  range  of  4.5  to  8 ym  of  the  approximate  distribution.^^ 
However,  this  same  condition  apparently  decreases  for  the 

larger  wavelengths. 

The  first  test  distribution  (Figure  5)  has  a mode  radius  that 
is  close  to  some  cloud  models;  however,  it  differs  in  the  concentration  of 
particles.  Comparing  the  results  for  the  extinction  coefficient 
with  Deirmendjian's  cloud  model, the  computed  values  were  slightly 
higher  (within  25*0)  for  the  comparable  wavelengths  of  0.45  ym,  1.61 
ym,  3.90  ym,  and  10  ym.  The  higher  values  are  to  be  expected  since 
the  mode  radius  of  Deirmendjian's  model  is  lower  (4  ym)  than  that 
of  the  test  distribution  (5  ym) . 

Comparing  the  results  to  the  "alto-stratus"  cloud  of 
Carrier,  et  at.,  which  is  closer  in  mode  radius  (4.5  ym) , the  com- 
puted values  of  are  lower  by  approximately  17  per  cent  for 

X = .488  ym,  1.06  ym  and  4.0  ym  and  30  per  cent  for  X = 10.6  ym, 
which  may  be  attributed  to  the  large  concentration  of  particles  in 
the  10-13  ym  range  in  the  "alto-stratus"  cloud.  The  computed  values 
for  8 (it)  are  again  low  by  42  per  cent,  5 per  cent,  54  per  cent,  and 
28  per  cent  for  X = .488  ym,  1.06  ym,  4.0  ym,  and  10.6  ym  respectively. 
With  the  exception  of  the  values  for  0.5  ym  and  4.0  ym,  these  lower 
values  may  similarly  be  attributed  to  the  broadened  concentration 

p.  1209. 

^®Deirmendj ian.  Scattering  and  Polarization,  p.  9. 
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above  10  vim  in  the  alto-stratus  distribution.  The  reason  for  the 
low  value  at  these  wavelengths  is  not  immediately  evident;  how- 
ever, Carrier  suggests  that  broadening  can  decrease  3(n)  by  factors 
in  the  SO  per  cent  range  and  so  these  values  are  not  considered  in- 
valid. I 


CHAPTER  IV 


ANALYSIS  OF  DATA 


The  data  produced  hy  the  computer  program  is  tabulated  in 
Appendix  1).  Tables  b and  7 contain  listings  of  the  extinction-  and 
backscattering-cocfficients  for  the  wavelengths,  step  size,  and  dis- 
tributions shown. 

To  determine  the  optinuuii  operating  conditions  we  consider 
the  wavelength  and  the  length  of  the  scattering  volume,  ignoring  for 
the  moment  the  practical  side  of  what  is  currently  available  in  the 
way  of  transmitter  power  and  special  requirements  ot  sensitive  de- 
tectors (such  as  low  temperature  operation) . Our  aim  is  to  maximize 
the  ratio  of  intensity  received  to  intensity  transmitted. 

The  wavelengths  utilized  were  selected  with  an  eye  towards 
the  transmission  "windows"  that  occur  in  the  infrared  region  of  tlic 
spectrum.  As  indicated  by  liquation  (1-b).  four  wavelength  dependent 
factors  affect  the  intensity  ratio:  (1)  the  field  of  view,  (2)  the 
extinction  factor  within  the  turbulence  itself,  (3)  the  :imount  of 
radiation  scattered  back  at  the  receiver,  and  (4)  the  transmissivity 
of  the  intervening  medium.  The  first  three  factors  have  been  in- 
cluded in  the  computer  program  computations  with  the  results  plotted 
for  individual  wavelengths  of  the  receiver  intensity,  normalized  to 
a transmitter  spectral  irradiance  of  1 watt/nf’A'crsus  the  distribution 
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Table  6.  Extinction  Coefficients  of  the  Various  Distributions, 


X(Ar) 

Distribution 

(micrometers) 

#1 

#2 

#3 

#4 

10.6  (.253) 

3.91808-02 

(0.170)* 

1.30680-01 

(0.568) 

4.73336-01 

(2.056) 

3.00532-01 

(1.305) 

10.6  (.506) 

3.91785-02 

(0.170) 

1.30681-01 

(0.568) 

4.73320-01 

(2.056) 

3.00427-01 

(1.305) 

4.0  (.255) 

8.46463-02 

(0.368) 

2.40473-01 

(1.044) 

4.32330-01 

(1.878) 

2.46862-01 

(1.072) 

4.0  (.477) 

8.47518-02 

(0.368) 

2.41875-01 

(1.050) 

4.34119-01 

(1.885) 

2.47657-01 

(1.075) 

1.06  (.253) 

6.80013-02 

(0.295) 

2.06574-01 

(0.897) 

3.99796-01 

(1.736) 

2.31788-01 

(1.007) 

0.50  (.251) 

6.53955-02 

(0.284) 

1.99431-01 

(0.866) 

3.90475-01 

(1.696) 

2.27353-01 

(0.987) 

•Figures  in  parentheses  indicate  loss  in  db/meter. 
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Table  7.  Backscatter  Coefficients  of  the  Various  Distributions, 
0(ir)  (meter) ‘^(steradian)"!. 


Distribution 


(micrometers) 

#1 

#2 

#3 

#4 

10.6 

(.253) 

6.40959-05 

1.99547-04 

2.1000-04 

1.00082-04 

10.6 

(.506) 

6.40957-05 

1.99548-04 

2.09992-04 

1.00032-04 

4.0 

(.255) 

1.24974-03 

3.87723-03 

1.18300-02 

3.37952-03 

4.0 

(.477) 

1.14929-03 

3.63343-03 

1.23361-02 

3.54947-03 

1.06 

(.253) 

3.59016-03 

1.05844-02 

2.31075-02 

1.21906-02 

0.50 

(.251) 

3.49286-03 

8.69341-03 

2.49455-02 

1.44447-02 

functions  in  Figures  6-9.  Figures  10-1.'.  are  comparison  plots  of 
the  normalized  receiver  intensity  for  the  different  wavelengths 
versus  a particular  distribution  function.  When  interpreting 
Figures  6-12,  it  is  important  to  remember  that  the  spectral  irradiance 
of  the  pulse  is  the  factor  that  has  been  kept  constant,  which  re- 
sults in  the  decrease  in  receiver  intensity  for  the  longer  lengths 
of  the  scattering  volume.  If  the  amplitude  of  the  pulse  is  kept 
constant  and  the  pulse  varied,  then  the  receiver  intensity  will  in- 
crease for  longer  length  of  the  scattering  volume.  The  transmis- 
sivity factors  for  all  four  wavelengths  are  computed  from  Table  2. 
From  the  calculations,  it  is  evident  that  the  property  of  high 
transmissivity  that  so  enhances  the  use  of  the  10.6  pm  laser  also 
hinders  it  in  this  study  as  the  backscattered  radiation  is  propor- 
tionately low.  From  Figures  10-12  it  is  clear  that  the  intensity 
ratio  of  the  4.0  pm  laser  is  significantly  greater  than  the  intens- 
ity ratios  of  either  the  10.6  pm,  the  1.06  pm,  or  the  0.5  pm  laser 
except  for  Distribution  #2,  which  indicates  that  the  intensity 
ratios  of  the  10.6  pm  line  is  comparable  to  that  of  the  4.0  pm  line, 
especially  when  the  transmissivity  of  the  two  are  taken  into  consi- 
deration. Furthermore,  the  extinction  of  the  0.5  pm  line  decreases 
its  intensity  level  to  a point  that  virtually  eliminates  it  from 
further  consideration. 

The  length  of  the  scattering  volume  is  important  since  it 
effectively  determines  the  amount  of  backscattered  radiation  pro- 
duced. Our  primary  concern  is  to  limit  the  backscattered  radiation 
to  that  which  is  singly  scattered  (thus  eliminating  multiply-scattered 
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velocities)  while  maximizing  the  amount  of  backscattered  intensity. 
Presumably,  an  unlimited  length  of  the  scattering  volume  would  re- 
sult in  the  maximum  return  of  unpolarized  radiation  from  all  modes 
of  scattering . To  determine  the  point  at  which  second-  and  higher- 
order  scattering  become  significant  we  relate  that  point  to  the 
optical  depth,  t = results  of  Anderson  and  Browell,^ 

it  appears  that  the  ratio  of  second-order  backscatter  to  first-order 
backscatter  is  about  10  per  cent  at  x = 0.25;  28  per  cent  at  a 

T of  1.0;  and  36  per  cent  at  a x of  4.0  for  a wavelength  of  0.9  vim. 
If  we  assume  similar  figures  for  the  wavelengths  of  1.06-,  4.0- , and 
10.6-ura  then  a reasonable  optical  depth  would  be  x = 1.0.  (We  note 
at  this  point  that  polarizing  our  initial  radiation  will  aid  in 
eliminating  at  least  50  per  cent  of  the  multiply-backscattered  radia- 
tion while  not  affecting  the  first-order  scattering.)  The  resultant 
length  of  scattering  volume  for  an  optical  depth  of  1.0  are  listed 
in  Table  8 for  each  of  the  wavelengths  and  distributions  considered. 

Distribution  #3,  corresponding  to  a concentration  of  500 
particles/cm  and  a mode  radius  of  10  vim,  is  seen  to  be  considerably 
lower  than  the  other  distributions  indicating  heavy  losses  for  this 
distribution  as  is  confirmed  by  Figures  6-9. 

Setting  the  depth  of  scattering  in  accordance  with  the  values 
in  Table  8 will  automatically  include  the  maxima  of  the  backscat- 
ter curves  in  Figures  5-10  since  the  maximum  (as  computed  by  setting 

^ Anderson  Browell,  First- and  Seaond-Order  Backsaatterina” 
p.  1348-1350. 
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Table  8.  Proposed  Single  Scattering  Lengths.* 


Wavelengths 

Distributions 

(micrometers) 

#2 

#5 

#4 

10.6 

7.64 

2.12 

3.32 

4.0 

4.16 

2.32 

4.04 

1.06 

4.88 

2.5 

4.32 

0.50 

5.01 

2.56 

4.39 

•Lengths  in  meters. 
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the  first  derivative  of  the  function  equal  to  zero)  of  the  Back- 
scatter  function  occurs  at  x = (2^^j.T^  which  is  % each  value  indi- 
cated in  Table  8.  We  relate  this  to  pulse  length  in  the  following 
manner;  a depth  of  scattering  of  four  meters  would  correspond  to 
a pulse  length  of  4m/3xl0®  rasec'^  = 1 I/S^IO"®  sec  or  about  thir- 
teen nanoseconds  which  has  been  achieved  in  lasers  of  wavelength 
2 

1.06  pm  and  4.0  pm,  and  can  theoretically  be  achieved  by  mode 
locking  or  modulating  by  absorbable  dyes  in  a 10.6  pm  laser. 

In  Chapter  II  we  computed  values  for  the  transmissivity  of 
the  atmosphere  (the  intervening  medium)  and  signal  power  for  the 
various  wavelengths.  Using  a metlian  value  of  transmissivity,  signal 
power  corresponding  to  an  SNR  = 10,  and  the  results  of  this  chapter, 
we  have  compiled  in  Table  9 representative  figures  for  the  required 
power  of  a laser  employing  a 6-inch  (15.24  cm)  Ions  (diffraction-limited 
detector),  operating  over  a 10  km  path  length.  The  figures  were  com- 
puted by  the  following  formula:  P.j.  = Req'd  Output  Powers  = 

AtPs 

-= (watts)  where  A.^  is  area  of  the  lens,  P^  is  the  re- 

T^(Scatter  Ratio) 

quired  signal  power,  is  the  transmissivity  of  the  atmosphere  and 
scatter  ratio  is  the  level  of  ackscattered  radiation  computed  for  a 
spectral  irradiance,  11(A),  of  I watt/meter^  over  the  scattering  lengths 

given  in  Table  8.  These  values  clearly  show  that  the  primary  wave- 
lenths  to  consider  are  10.6  pm  and  1.06  pm  with  the  former  having  a 
slight  edge  in  average  power  required. 

2 A.  Kestenbaum.  "Laser  Processing  Fundamentals."  Lasers  in 
Industry,  ed.  S.  S.  Charschan.  New  York:  Van  Nostrand  Reinhold 
Co.,  1972. 
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Table  9.  Power  Required  for  Operation  of  the  System* 


Wavelength 

(micrometers) 

10.6 

4.0 

1.06 

0.50 

Transmissivity  (t^) 
(dimensionless) 

4.7x10"5 

6.2  xl0"8 

4.4  xlO'3 

8.9  xl0"5 

Required  Signal 
Power  (watts) 

1.41xl0'll 

9.94x10"^! 

1.41x10"^ 

6.36x10‘3 

Scatter  Ratio 
(dimensionless) 

Distribution  #2 

1.30x10"^2 

1.85x10"^2 

4.05xl0'13 

8.1  xl0"l** 

Distribution  #3 

3.6  xl0‘13 

3.5  xl0‘12 

ro 

O 

X 

00 

1.18x10‘13 

Distribution  #4 

2.8  xl0'13 

1.6  xl0"12 

4.4  xl0"13 

1.18x10"^3 

Required  Output 
Power  (watts) 

Distribution  #2 

4.21x103 

1.58x102 

1.44x10** 

1.61x10® 

Distribution  #3 

1.52x10'* 

8.36x10^ 

1.22x10** 

1.1  xlo2 

Distribution  #4 

1.95x10'* 

1.83x102 

1.33x10** 

1.1  xlo2 

Avei  ge 

1.3  xio** 

1.42x102 

1.33x10** 

7.87x10^ 

‘Average  power. 
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Normalized  releiver  intensity,  watts 


Normalized  receiver  intensity,  watts 


Nornalized  receiver  intensity,  watts 


Normalized  receiver  intensity,  watts 
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3 

Length  of  Scattering  Volume,  Meters 


Figure  10.  Plot  of  normalized  receiver  intensity  vs.  length  of  scattering 
volume  at  all  wavelengths  for  distribution  #3. 


Length  of  Scattering  Volume,  Meters 


Figure  11.  Plot  of  normalized  receiver  intensity  V5.  length  of  scattering 
volume  at  all  wavelengths  for  distribution 


Normalized  receiver  intensity,  watts 


CHAPTER  V 

SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATION 

The  purpose  of  this  paper  was  to  determine  the  feasibility 
of  utilizing  pulsed  lasers  to  determine  the  velocity  structure  of 
the  turbulence  associated  with  tornadoes.  To  this  end,  a computer 
program  was  developed  using  Mie  parameters  to  obtain  the  extinction 
and  backscatter  properties  associated  with  the  particles  presumed 
to  bo  distributed  within  the  tornado.  The  program  was  then  used  to 
analyze  three  possible  distributions  of  particles  with  four  different 
wavelengths.  The  data  obtained  were  (1)  tlie  extinction  and  backscatter 
coefficients  of  the  various  wavelengths  for  those  distributions,  and 
(2)  the  associated  backscatter  intensity  ratios  (for  the  various 
distributions  and  wavelengths)  evaluated  per  length  of  scattering 
volume  for  a specified,  diffraction  limited  f ield-of-view  and  a 
specified  scattering  area.  In  addition,  the  computer  program  will 
output  the  Mie  parameters  for  any  desired  value  of  the  size  parameter, 
X ■ 2itr/X,  and  it  is  adaptable  to  other  continuous  distributions  or 
a combination  of  continuous  distributions. 

It  was  determined  that  a SNR  = 10  could  be  achieved  using 
pulse  widths  of  the  order  of  10-20  nanoseconds  and  output  powers  of 


the  order  of  10  kilowatts.  To  achieve  these  figures,  we  used  a 10 
km,  clear-air  approach  medium  and  an  aligned,  diffraction-limited. 
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heterodyne  receiver  with  a 15.24  cm  (diameter)  lens. 

Of  the  wavelengths  considered,  it  was  determined  that 
10.6  ym  and  1.06  ym  are  the  most  conducive  to  this  type  of  use 
while  0.5  ym  will  be  severely  handicapped  by  the  funnel  cloud 
distributions  and  4.0  vm  is  a victim  of  the  intervening  medium 
(specifically,  molecular  absorption). 

Three  distributions  (designed  as  limiting  distributions), 
with  particle  radii  on  the  order  of  1-25  ym,  were  utilized  to  re- 
present probable  conditions  found  in  tornado  funnel  clouds.  The 
results  obtained  were,  in  general,  applicable  to  all  three  distri- 
butions . 

Conclusions 

From  the  data  obtained  in  this  study,  it  appears  that  a 
pulsed  laser  operating  in  the  near  infrared  has  the  capability  to 
deliver  sufficient  power  to  project  electromagnetic  radiation 
through  the  atmosphere  to  a range  of  up  to  10  km,  scatter  that  radia- 
tion from  a turbulent  distribution  of  particles,  to  return  and  to 
provide  a measureable  signal  to  a coaxially  mounted  detector.  Table 
10  indicates  the  power  required  for  various  ranges  up  to  10  km. 

The  velocity  spectra  obtained  as  a result  of  implementing 
such  a system  is,  of  itself,  limited  in  the  sense  that  a velocity 
so  obtained  is  the  total  instantaneous  velocity  of  a particle  in 
the  direction  of  the  observer  at  the  time  of  scattering.  There  is 
no  discrimination  between  the  various  translational  and  rotational 
components  of  this  velocity.  Further,  there  is  no  indication  of 
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Table  10.  Minimum  Power  Requirements  at  10.6  pm  for  the  Various 
Distributions  versus  Range  (Power  in  Watts).* 


Range  (kilometers) 


Distribution 

1 

3 

5 

8 

10 

#2 

0.54 

3.93 

28.9 

573.4 

4.2x103 

#3 

1.94 

14.2 

104.2 

2.07x103 

15.2x103 

#4 

2.49 

18.3 

134.0 

2.66x103 

19.5x103 

Average 

1.66 

12.1 

89.0 

1.77x1o3 

13.0x103 
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movement  in  the  direction  vertically  perpendicular  to  the  direction 
of  propagation.  However,  it  is  anticipated  that  one  may  obtain 
some  of  this  information  by  moving  the  laser  beam  in  a coordinated 
manner  both  horizontally  and  vertically  a,id  thus  obtain  a mapping 
of  the  instantaneous  velocities  of  the  tornado.  Comparison  between 
the  velocities  obtained  on  one  side  of  the  tornado  and  its  horizon- 
tal counterpart  on  the  other  side  of  the  tornado  would  provide  some 
indication  of  the  relationship  between  the  rotational  and  transla- 
tional velocities  since  the  two  spectra  then  represent  appro..iching 
and  receding  velocities. 

The  data  obtained,  though  limited  to  the  distributions 
cited,  is  considered  a good  indication  of  the  limiting  conditions 
likely  to  be  found  in  a funnel  cloud.  There  are  other  conditions 
that  are  natural  extensions  of  the  work  initiated  here  and  resulting 
analysis  of  these  conditions  may  be  conducted  utilizing  the  existing 
computer  program.  Examples  include  the  investigation  of  distribu- 
tions that  would  result  from  contact  with  the  ground  such  as  the  dis- 
tribution of  prairie  dust  or  that  of  a plowed  field.  Also,  investi- 
gation of  the  distribution  of  rain  should  be  analyzed  to  determine 
to  what  extent  the  signal  is  further  degraded.  Probable  distribu- 
tions of  rain  have  been  suggested  by  the  work  of  Rensch  and  Long.^ 

It  is  anticipated  that  rain  could  conceivably  eliminate  the  possi- 
bility of  sufficient  backscatter  being  received  by  the  detector  and 

^Rensch  and  Long,  Comparative  Studies,  pp.  1566-67. 
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the  need  to  select  an  optimum  direction  from  which  one  would  direct 
the  laser  beam  becomes  apparent. 
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Recommendation 


A.1  . ' ■' 

Laser  technology  has  progressed  to  the  point  that  velocity 


discrimination  on  the  order  of  1 m/sec  is  achievable  outside  the 
laboratory.  Furthermore,  detectors  are  available  (at  the  cited 
wavelengths)  with  sensitivities  of  the  order  required  to  process 
the  returned  signal.  However,  it  is  recommended  that  a system  con- 
taining the  required  components  including  the  ancillary  equipment 
be  constructed  and  tested  to  verify  the  results  obtained  in  this 
study . 


APPENDIX  A 


ADAPTATION  OF  MIE  PARAMETERS  TO  COMPUTER  ITERATIONS 


The  objective  is  to  obtain  the  expressions  (2-5)  through 
(2-8)  in  the  text.  Therefore  we  must  express  a^,  b^, 


and  A in  iterative  form  and  determine  the  initial  values  of 
n 


those  iterations.  These  are, 
A 


(-2-+-)R  (W  }-R  OV  ,} 
*■  m X e n e ' n-1 


n 


(-2-+-)W  -W  , 
' m X n n-1 


(A^m.5)Re{»V^}-R^{W^.l} 


^ (A  m + -)W  -W  , 
'■  n x"^  n n-I 


W = w - W - 

n x n-1  n-2 


W = sinx  + icosx 
o 


W 2 = cosx  - isinx 


, _ n 1 

n ~ ” mx  * n 


— - A , 
mx  n-1 


. ^ cos  mx 

Ao  - = liiTisr 


Since  mx  is  a complex  function,  we  expressed  this  value  in  the 
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following  form:  with  m = r + ik 

. _ cos  rx  sin  rx  - i sinh  kx  cosh  kx 

Ap  a 

sin^rx  cosh^kx  + cos^rx  sinh^kx 

instead  of  the  value  used  by  Deimendjian  in  the  cited  reference. 

The  values  for  the  angle  dependent  function  and  were 
computed  for  0 = ir  prior  to  inserting  them  into  the  program.  Their 
values  are: 


i 


^Deirmendjian  and  Clasen,  Light  Saattenng^  pp.  35-36. 


APPENDIX  B.  FLOWCHART  OF  COMPUTER  PROGRAM 


Note 


Compute 


ext*  sc* 


r 
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Note  1 


Note  2 


Note  3 


Note  4 


Note  5 


71 


Number  of  array  elements,  n,  is  determined  by  size  para- 
meter X.  As  a minimum,  n = 1.2x  + 12. 

Wavelength  is  expressed  in  micrometers;  Index  of  refrac- 
tion is  a complex  function.  The  values  used  in  this  paper 
are  listed  in  Table  Bl. 

The  Mie  Coefficients  and  Iteration  variables  are  computed 
for  each  x as  described  in  Appendix  A. 

The  Mie  Parameters  are  computed  for  each  x as  described  in 
the  text  and  repeated  for  convenience  below.  Note,  how- 
ever, 

2 1.2X+12 

K fx,m)  = — 7 7 (2n+l)Re{a  +b  } (dimensionless) 
x^  n n 


1.2X+12 

' :::2  I r*|b„| 

^ n = l 


(dimensionless) 


1.2X+12 


Sj(x,m,n)  = -S2(x,m,Tr)  = [ 

n=l 


(dimensionless) 

that  is  not  required  for  our  study  and  was  not  speci- 
fically mentioned  in  the  text.  It  is  the  normalized  Mie 
scattering  cross  section  and  is  an  added  feature  of  the 
program.  One  may  wish  to  use  it  for  various  other  func- 
tions including  the  calculation  of  the  normalized  Mie  ab- 
sorption cross  section. 

The  distribution  functions  are  of  the  form 


t 
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fiU.A)  = a(^)^*Vexp{-b(^)Yx7} 

j with  the  constants  a,  b,  a,  y given  in  Table  4 , for  the 

various  distributions,  i. 

' Note  6:  The  extinction-,  scattering-,  and  backscatter-coefficients 

are  computed  for  each  distribution  by  numerical  integra- 

Ition.  The  formulas  are  given  below  for  convenience.  See 

Note  4 in  regards  to  the  scattering  coefficient. 

^2 

ax^^'^expf-bC^)  V}K.g^^dx  m"! 

^1 

®sc^“’^’*l’*2^  = ^"^“(10-6)1  ^ax“'^^exp{-b(^)\''^}Kg^dx  m-l 

*2 

S(ir,X,m,Xj,X2)  = j ax“exp{-b(^)^xY}  |Sj(v)  |^dx  m"l-sr- 1 

Note  7:  The  intensity  ratios  (BACKSCATTER)  are  computed  for  a 

specific  field  of  view  and  scattering  area  (OMEGA)  over 
values  of  the  scattering  length  (RESOLUTION)  of  from  1 
to  10  meters  according  to  the  formula. 

RECEIVER^I^NTENSITY  ^ AL0B(TT)exp(-2LBg^^)  (dimensionless) 

Note  that  it  does  not  consider  the  transmissivity  of  the  intervening 
medium. 


I : 

j 


1 


Table  Bl.  Indices  of  Refraction 


Wavelength  Complex  Index  of 

(Micrometers)  Refraction,  RM(=m) 


10.60 

1.212-i 

0.0601 

4.00 

1.353-i 

0.006 

1.06 

1.33-1 

0.000 

0.50 


1.3A-i  0.000 
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31 

32 

33 

34 

35 

36 
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51 
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54 
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W( 2)=CHPLX( SXN.CXN) 

SlNPaSINlP ) 

COSP*COS(P) 

SlNHQ=(eXP(Q)-EXP(-Q) »/2 
C0SHQ3(eXP<Q)>-EXP(-Q)  )/2 

01 V=( ( SI NP**2 ) *(C0SHQ**2  > » ♦( i COSP * *2 ) * ( S 1 NMQ* *2 » I 
RMR=CQSP*S INP/DI V 
RMls>( SINHOOCOSHO/OIVI 
0<  2»=CMPLX<RMR,RMI  1 
DO  15C  N=3,M,1 

W(  N)  = ( ( ( 2*  ( N-2  ) )-l  ) ♦V*(N-l  ) /XN)  -WIN-a  » 

D<N»  = < - { N-  2)/RMX)<-  ( I/(  ( ( N-2  ) /RMX)-0(N>1  ) I ) 
RWN=REAL(W(N) ) 

RWNl  =REAL(  W(  N-n  ) 

A(N»  = ( RWN*  U (N-2) /XN  )4-(0(N  )/RH  ))  -RWNl  ) / 

*<W<N)  ♦(  ( (N-2  )/XN)«-(D  (N  )/RM  ) ) - W(N-1  ) ) 
a(  N)  = ( RWN*  ((  (N-2  ) / XN)  4'(0(N  )*RM  ) )-RWNl  )/ 

♦ (W(N) *( ( (N-2  )/XN ) + (D{N  )*KM  ) )-W (N-1  ) ) 

PI=-((-l  )**(N-2)  )★(  (N-2)*(N-1  ) )/2 

AMP=(  (2.4(N-2))^1)4>P14(A(N)-8(N))/((N-2)*(N-1)) 
RAMP=RAMP* AMP 
ANBN  = A(  N)4-B(  N) 

RAN8N=REAL( AN8N) 

SKEXTN=(  2*  (N-2  )-»-l  ) *RAN8N 
SKEXT  = SKEXT<-SKEXTN 
RsCA8S( A(N ) ) 

S=CABS(6(N) ) 

SKSCN=(2*(  N-2  )*-l  )4(R**2*S4*2) 

SKSC=SKSC'*-SKSCN 
150  CONTINUE 

RAMPN=CABS (RAMP) 

RINT^O . 1 5* (RAMPNO42.0 ) 

EXT=2.C*SKEXT/ (XN**2 .0 ) 

SCT=2,C*SKSC/(  XNW*2.0) 

R=XN*U/( 2.0*3. 1415926) 

XOKl=R** 1. 33333 
XOK23=R**2 .0 
XOK4sR**3.0 
XDK5=R**5.0 
XDK6sR*«6.0 
XOKe=R*«a.O 
AA*0.1 1098E-06 
BB*0.7901E-03 

FUNl  = 2 3.4a02  32*XOK4*EXP( -0 .423  41 *X OK  1 ) 

FUN2*1 .0*XOK5*EXP( -0 . 1 *XOK23) 

FUN3  30 .C  0781 24  8*X0K5*EXP( -0.02  5*XOK2  3) 
F0N4=AA*XDK8*EXP( -88*XDK4) 

HI «( (U/( 2.0*3. 1415926) )**3.3 )*(1 .C0E-C6) 

H23H1 
H4  «H1 
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1 

CONI  32.0 *3. 141 5926*0 .154SKEXT 

1 80 

CON232.0*3.141S926«C .IS«SKSC 

1 ®* 

FIRSr«XN-0 .01 

62 

LASTl sXN-9.61 

83 

LAST2=XN-I4.7l 

84 

IF  (FIRST)  200.160.155 

as 

155 

IFCLASTl ) 180.1 oS. 158 

86 

158 

IF(LAST2 1182,170.230 

67 

160 

FUNl=FUNl/2,0 

88 

FUN2=FUN2/2.0 

89 

FUN3=FUN3/2.0 

90 

FUN4=FUN4/2. 0 

91 

GO  TO  180 

92 

165 

FUNl=FUNl/2.0 

93 

GO  TO  180 

94 

170 

FUN2=FUN2/2.0 

95 

FUN3=FUN3/2,0 

96 

FUN4=FUN 4/2.0 

97 

GO  TO  182 

98 

180 

BEXT1=BEXT1*(C0N1*FUN1*HI ) 

99 

BSCl=dSCl+(CON2*FUNl *H1 ) 

100 

BPIl =BPI I* (RINT*FUN1 ♦HA) 

101 

182 

BEXT2=BEXT2*(C0N1 ♦FUN2*H2) 

102 

BSC2*OSC24(CON2*FUN2*H2) 

103 

8PI2  = t3PI2«-  (RINT*FUN2*H4  ) 

104 

BEXT3=BEXT  3* ( C ON  1 *FUN3*H2 ) 

lOS 

BSC33BSC3* (C0N23FUN3  *H2 ) 

106 

BPI3=8PI3*(RINT*FUN3*H4) 

107 

BEXT4=BEXT4*(C0N1 ♦FUN4*H2) 

\ 106 

BSC4=BSC4* (C0N2*FUN4*H2) 

1 109 

BPI43BPI4* (RINT*FUN44H4) 

110 

00  195  N3| 51 . 1 50 1 . 150 

111 

YN=N*0.C 1 

112 

DIFF3XN-YN 

113 

IF<OIFF)200,185,195 

114 

185 

WR  ITE( b, 190) XN.eXT ,SCT .RAMP 

IIS 

190 

F0RMAT(5E14,6//) 

116 

GO  TO  200 

117 

195 

CONTINUE 

1 lie 

20  0 

CONT INUE 

119 

00  300  J3l  .4 . 1 

12'J 

«RITE(6.205) 

^ 121 

20  5 

F0RMAT(4X.«J»,5X.»BEXT«,10X,*BSC».l 

1 X,  'BIPI ) •/ ) 

;i  122 

GO  TO  (2 10 .220  .230  .240  ). J 

123 

21  0 

WRITE(6,21 5) J.dEXTl.BSCl .BPI 1 

124 

21  5 

FORMAT ( 2 X.  I 3 , 2 X . 3E  1 4 .6/ ) 

125 

BEXT=3EXTl 

126 

BETA3BP I 1 

127 

GO  TO  280 

il 
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128 

220 

WR1TE(  6*225) J.BEXT2.BSC2.DP12 

129 

22  5 

FORMAT(2X«  1 3 *2 X . 3E 14 .6/ ) 

130 

BEXT=BEXT2 

131 

BETA=BPI 2 

132 

GO  TO  280 

133 

230 

W«ITE(6.235) JtDEXT3,BSC3.BPl3 

134 

235 

FORMAT(2X.  I3.2X,3El4.b/) 

1 35 

BEXT=BEXT3 

136 

BETA-BPI 3 

137 

GO  TO  280 

136 

240 

WRITE(6.245) J.BEXT4,8SC4,BPI4 

1 39 

245 

FORMAT {2X, I 3 , 2 X , 3E 1 4 .6/ ) 

140 

BEXT=8EXT4 

141 

BETA-BPM 

142 

260 

WRITeC6.282) 

143 

262 

FORMAT {/3X , • RESOLUTION* .4X , • SOLID  ANGLE* . 3X, • BACKSCATTER 

144 

DO  290  L=l .10,1 

145 

RESO=L 

146 

OMEGA=0.616E-Oa 

147 

BACKSC=RESa*0MEGA4BETA*EXP{-2*RES0*>BeXT) 

146 

•RITE  (6,285)RES0,0MEGA,BACK3C 

149 

285 

FORMAT(/3E 14.6/) 

150 

OELTA=BACXSC-l .OOE-SO 

151 

IF (DELTA) 300 .300.290 

152 

290 

CONTINUE 

153 

30  0 

CONTINUE 

154 

STOP 

155 

END 

SEXEC 

X 

KEXT  KSC  SI 

0.151000E  01 

0.3ai408E  00  0.129208E  00  0.8265195-01  0.116796E  CO 

0.301000E 

01 

0. lCa916E 

01 

0. 

597752E 

00  -0. 

185405E 

00 

0. 159753E-ai 

0.4S1000E 

01 

0. 181213E 

01 

0. 

1 1 4461E 

0 

1 

O 

• 

10422CE 

00 

a.44902SE-01 

0.601000E 

01 

0. 237685E 

01 

0. 

157640E 

01  -0. 

1417  1 IE 

00 

-0.4381S8E  CO 

0.7S1000E 

01 

0.270629E 

01 

0. 

181043E 

01 

0.304792E 

CO 

0. 16e074E 

00 

0.901000E 

01 

0.281  91  IE 

01 

0. 

184958E 

01 

-0. 253443E 

00 

-0.307156E 

00 

79 


0.1051006 

02 

0.2767336 

01 

0.1744176 

01 

0. 6846216 

00 

0.273401 6 

0.1201006 

02 

0. 2626056 

01 

0.1564 766 

01 

-0.5238916 

00 

0.5193016 

0.13S100E 

02 

0.2466906 

01 

0.1379226 

01 

0.6043956 

00 

0.2123896 

J BEXT  8SC  8(PI) 

1 0.3918086-01  0.2294746-01  0.6409596-04 


RESOLUT ION 

SOLID  ANGLE 

BACKSCATTE9 

0.1 OCCOCE 

01 

0.6160006-08 

0.365C72E-12 

0.2000006 

01 

0.6160006-08 

0.6751 136-12 

0.3000006 

01 

0.6160006-08 

0.9363456-12 

0.4000006 

01 

0.6160006-08 

0.11 54  366-  1 1 

0.5000006 

01 

0.6160006-08 

0.1334206-1  1 

0.6000006 

01 

0 .6160006-08 

0.  1480376- 1 1 

0.7000006 

01 

0.6160006-08 

3. 1596926-1 1 

0.8000006 

01 

0 .6160006-08 

0.1687506-  1 1 

0.9000006 

01 

0.61b00CE-08 

0.1755356-11 

0.10C00C6 

02 

0.6160006-08 

0.1803396-11 

J BEXT 

SSC 

B(PI  ) 

2 0.1306806  00  0.7774076-01  0.1995476-03 


4 


ik 


RESOLUT ION 

SOLID  ANGLE 

BACKSCATTER 

OalOCOOCE 

01 

0.616000E*C8 

0.946496E-12 

0*200000E 

01 

0.616000E>08 

0. 14S761E-1 1 

0.300000E 

01 

C.616000E-08 

0. 168355E-1 1 

0.400000E 

01 

0.616000E-08 

0.172e45E-l  1 

0.500000E 

01 

0. 616000 E>08 

0. 166364E-1  1 

0.60C000E 

01 

0.6160CCE>0a 

0.153721E-11 

0.70000CE 

01 

0 .OlBOOOE-OO 

0.138C93E-11 

0.60000CE 

01 

0.61 6000E- 08 

0.121S23E-11 

0.900000E 

01 

0.6160COE>08 

J. 10  5269e-l  1 

O.IOOOOOE  C2  C.61600CE-08  0.900642E-12 

J BEXT  BSC  a(PI) 

3 0.473336E  00  3.31C62IE  00  0.210000E 

RESOLUTION 

SOLID  ANGLE 

BACKSCATTER 

O.IOOOOOE 

01 

0.6160C0E-0a 

0.S01957E-12 

0.200000E 

01 

0.61600CE>08 

0.3895S0E-12 

0.3000COE 

01 

0.616000E-08 

0.226736E-12 

!i 


61 


0«400000E 

01 

0.616000E>0a 

0.11 7307E-12 

0.50000CE 

01 

0.616000E-0a 

0.S68984E-I3 

0.60000CE 

01 

0.6I60COE-Oa 

0.264940E- 1 3 

0.700000E 

01 

0.61600CE-0ei 

0.11 99 J9b- 1 3 

OcSOOOOCE 

01 

0 .616000E-0a 

O.S31Bb5E>14 

e.900000E 

01 

0 .oiaocoE-os 

0.2321 Qot-l * 

O.tOOOCCE  02 

j acxT 

4 0.3005 

0.e»1600CE-CH 

asc 

52E  CO  0.I92SM 

0. 1001Coe>  14 

81  Ml  1 

9E  00  0 . 1 00  3<l2r -C  3 

RESOLUT ION 

SOCIO  ANGLE 

BACKSCAT  TER 

0. 1 OCOOCE 

01 

0 .61 600CE-0a 

0.3J79B5E-12 

0.200000E 

Cl 

0.61 600CC- OS 

0. 370585k- 12 

0.300000E 

01 

0.61 6000E-0B 

0.30474BE- 12 

0.4000COE 

01 

0.616000E-0b 

0.2227628- 12 

j 

O.SOOOOOE 

01 

0.616000E-08 

3. IS26S5E- 12 

0.400C00E 

0 1 

0.616000E-0a 

0. 100428E>I2 

! 

0.700000E 

01 

0 .61600CE-0a 

0.642335E- 1 3 

i 

1 

o.aoooooE 

01 

0.616000E-08 

0.402453E-13 

0.900000E 

01 

0.616000E*08 

0.248215E-13 

1 

O.IOOOCOE 

02 

0.61600CE- 08 

0. 1S1198E>13 

'i 

i 

1 

! 1 

J 

APPENDIX  D 


The  tables  that  follow  are  the  output  of  the  computer  pro- 
gram. The  Mie  parameters  are  not  included  here  since  their  tabula- 
tion is  not  directly  related  to  the  major  emphasis  of  this  study. 

An  example  has  been  included  in  the  actual  computer  printout  con- 
tained in  Appendix  C. 

For  all  tables,  an  explanation  of  t«  .s  is  given  as 
follows:  Exponents  are  expressed  as  follows:  1.5*10*^  • l.S-0.2 
Resolution  ? length  of  scattering  volume 
Solid  Angle  i scattering  area  < field  of  view 


Rackscatter  i 


■receiver 


It  does  not  include  trans- 


iissivity  of  intervening  medium. 


A 


r 

\ 

t 


83 


Table  Dl.  Intensity  Ratios  for  Distribution  1;  Wavelength 
■ 10.6  im,  Ax  = 0.15;  Sgxt  “ 3.91808-02in'^; 

6^^  » 2.29474-02m'l;  e(n)  = 6.40959-05  m"l-sr'^ 


Resolution 

(Meters) 

Solid  Angle 
(Meter) 2 (Steradian) 

BACKSCATTER 

(Dimensionless) 

1.0 

6.16-09 

3.65072-13 

2.0 

tl 

6.76113-13 

3.0 

II 

9.36345-13 

4.0 

tl 

1.15436-12 

5.0 

II 

1.33420-12 

6.0 

If 

1.48037-12 

7.0 

II 

1.59692-12 

8.0 

II 

1.68750-12 

9.0 

II 

1.75535-12 

10.0 

If 

1.80339-12 

I 

i 


Table  D2.  Intensity  Ratios  for  Distribution  2;  Wavelength 
= 10.6  um.  Ax  = C.15;  = 1.30680-01 

= 7.77407-02  m-l;  B(ir)  = 1.99547-04  m'l-sr"l. 


Resolution 

(Meters) 

Solid  Angle 
(Meter) 2 (Steradian) 

BACKSCATTER 

(Dimensionless) 

1.0 

6.16-09 

9.46496-13 

2.0 

II 

1.45761-12 

3.0 

11 

1.68355-12 

4.0 

II 

1.72845-12 

5.0 

II 

1.66364-12 

6.0 

II 

1.53721-12 

7.0 

II 

1.38093-12 

8.0 

It 

1.21523-12 

9.0 

II 

1.05269-12 

10.0 


9.00642-13 


8S 


Table  D3.  Intensity  Ratios  for  Distribution  3;  Wavelength 

» 10.6  pm.  Ax  = 0.15;  6 ext  = 4.73336-01  m"!; 

= 3.10621-01  6(ir)  = 2.10000-04  m‘^-sr"^ 

sc  ’ \ j 


Resolution 

(Meters) 


Solid  Angle 
(Meter) 2(Steradian) 


6.16-09 


BACKSCATTER 

(Dimensionless) 


5.01957-13 


3.89550-13 


2.26736-13 


1.17307-13 


5.68984-14 


2.64940-14 


1.19939-14 


5.31885-15 


2.32186-15 


1.00106-15 
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Table  D4.  Intensity  Ratios  for  Distribution  4;  Wavelength 
= 10.6  ym.  Ax  = 0.15;  6ext  ~ 3.00532-01 

8^^  = 1.92589-01  m'l;  6(Tr)  = 1.00082-04  ra"l-sr‘l 


Resolution 

(Meters) 

Solid  Angle 
(Meter) ^ (Steradian) 

BACKSCATTER 

(Dimensionless) 

1.0 

6.16-09 

3.37895-13 

2.0 

It 

3.70585-13 

3.0 

II 

3.04748-13 

4.0 

It 

2.22762-13 

5.0 

M 

1.52655-13 

6.0 

tl 

1.00428-13 

7.0 

11 

6.42335-14 

8.0 

II 

4.02453-14 

9.0 

II 

2.48215-14 

10.0 

II 

1.51198-14 

Table  D5.  Intensity  Ratios  for  Distribution  1;  Wavelength 
= 4.0  um.  Ax  = 0.40;  Sgxt  “ 8.46463-02 
8 = 7.78943-02  8(11)  = 1.24974-03  m-^-sr"! 

SC 


Resolution 

(Meters) 

Solid  Angle 
(Meter) ^(Steradian) 

MCKSCATTER 

(Dimensionless) 

1.0 

8.77-10 

9.25329-13 

2.0 

tl 

1.56244-12 

3.0 

II 

1.97866-12 

4.0 

II 

2.22735-12 

5.0 

II 

2.35058-12 

6.0 

II 

2.38141-12 

7.0 

II 

2.34562-12 

8.0 

II 

2.26322-12 

9.0 

II 

2.14959-12 

10.0 

II 

2.01647-12 
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Table  D6.  Intensity  Ratios  for  Distribution  2;  Wavelength 
» 4.0  im.  Ax  = 0.40;  ® 2.40473-01 

8 * 2.18361-01  6(ir)  = 3.87723-03  m-^-sr"^ 

SC 


Resolution 

Solid  Angle 

BACKSCATTER 

(Meters) 

(Meter)2(Steradian) 

(Dimensionless) 

1.0 

8,77-10 

2.10208-12 

2.0 

ft 

2.59900-12 

3.0 

ft 

2.41005-12 

4.0 

It 

1.98652-12 

5.0 

It 

1.53508-12 

6.0 

If 

1.13878-12 

7.0 

It 

8.21324-13 

8.0 

If 

5.80276-13 

9.0 

11 

4.03567-13 

10.0 


2.77205-13 
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r 


[ 

I 

1 

Table  D7.  Intensity  Ratios  for  Distribution  3;  Wavelength 
■ 4.0  wm.  Ax  = 0.40;  8^^^  = 4.32330-01 

1 8 * 3.59059-01  8(11)  “ 1.18300-02  m“l-sr‘l 

1 sc 


Resolution 

(Meters) 

Solid  Angle 
(Meter) 2 (Steradian) 

BACKSCATTER 

(Dimensionless) 

1.0 

8.77-10 

4.36984-12 

2.0 

ft 

3.68111-12 

3.0 

II 

2.32569-12 

4.0 

It 

1.30609-12 

5.0 

If 

6.87648-13 

6.0 

II 

3.47560-13 

7.0 

11 

1.70789-13 

8.0 

II 

8.22118-14 

9.0 

If 

3.89556-14 

10.0 


1.82310-14 


I 


Table  D8.  Intensity  Ratios  for  Distribution  4;  Wavelength 

■ 4.0  pm.  Ax  * 0.40;  * 2.46862-01  m"'; 

8.^  - 1.96040-01  m-»;  6(n)  * 3.37952-03  m‘»-sr-l 

sc 


Resolution 

(Meters) 

Solid  Angle 
(Meter) 2 (Steradian) 

BACKSCATTER 

(Dimensionless) 

1.0 

8.77-10 

1.80897-12 

2.0 

fl 

2.20821-12 

3.0 

II 

2.02167-12 

4.0 

II 

1.64523-12 

5.0 

II 

1.25520-12 

6.0 

fl 

9.19334-13 

7.0 

II 

6.54633-13 

8.0 

II 

4.56634-13 

9.0 

II 

3.13544-13 

10.0 

It 

2.12635-13 
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Table  D9.  Intensity  Ratios  for  Distribution  1;  Wavelength 
■ 1.06  ym.  Ax  » 1.5;  ■ 6.80013-02 

B.  ■ 6.80013-02  6(it)  * 3.59016-03  m-I-sr-‘ 

SC 


Resolution 

(Meters) 

Solid  Angle 
(Meter) ^ (Steradian) 

BACKSCATTER 

(Dimensionless) 

1.0 

6.16-11 

1.93032-13 

2.0 

11 

3.36972-13 

3.0 

II 

4.41184-13 

4.0 

II 

5.13444-13 

5.0 

II 

5.60193-13 

6.0 

II 

5.86751-13 

7.0 

II 

5.97497-13 

8.0 

II 

5.96023-13 

9.0 

II 

5.85261-13 

10.0 


5.67600-13 
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Table  DIO.  Intensity  Ratios  for  Distribution  2;  Wavelength 

■ 1.06  win.  Ax  = 1.5;  B - 2.06574-01  m’h 

ext 

B - 2.06574-01  e(iT)  = 1.05844-02  m-l-sr'‘ 

sc 


Resolution 

(Meters) 

Solid  Angle 
(Meter) 2(Steradian) 

BACKSCATTER 

(Diaensionless) 

1.0 

6.16-11 

4.31341-13 

2.0 

II 

5.70719-13 

3.0 

II 

S. 66351-13 

4.0 

II 

4.99570-13 

5.0 

II 

4.13122-13 

6.0 

11 

3.27968-13 

7.0 

II 

2.53134-13 

8.0 

It 

1.91388-13 

9.0 

It 

1.42443-13 

10.0 


1.04705-13 


J 
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Table  Dll.  Intensity  Ratios  for  Distribution  3;  Wavelength 

■ 1.06  ym.  Ax  = 1.5;  ” 3.99796-01 

B = 3.99796-01  6(11)  = 2.31075-02  m'^-sr”^ 

sc 


Resolution  Solid  Angle  BACKSCATTER 

(Meters)  (Meter) 2 (Steradian)  (Oinensionless) 


1.0 

2.0 

3.0 

4.0 


5.0 

6.0 


7.0 

8.0 

9.0 


6.16-11 

It 


6.39845-13 

5.75237-13 

3.87864-13 

2.32466-13 

1.30621-13 

7.04586-14 

3.69507-14 

1.89826-14 

9.59955-15 


j 


1 


10.0 


4.79457-15 
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Table  D12.  Intensity  Ratios  for  Distribution  4;  Wavelength 
■ 1.06  pm.  Ax  = 1.5;  ® 2.31788-01  m"^; 

0^^  « 2.31788-01  B(ir)  = 1.21906-02  m"l-sr‘^ 


Resolution 

Solid  Angle 

BACKSCATTER 

(, 

1 

(Meters) 

(Meter) 2 (Steradian) 

(Dimensionless) 

1 

1.0 

6.16-11 

4.72363-13 

i. 

1 

i 

i: 

2.0 

ft 

5.94261-13 

r 

r 

3.0 

11 

5.60713-13 

i 

i. 

4.0 

It 

4.70274-13 

i 

1 

5.0 

tt 

3.69771-13 

i 

6.0 

tt 

2.79117-13 

1 

7.0 

tt 

2.04835-13 

i 

1 

8.0 

tt 

1.47254-13 

1 

I 

9.0 

It 

1.04206-13 

5 

10.0 


7.28318-14 


1 

1 


[ 

i 
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Table  D13.  Intensity  Ratios  for  Distribution  1;  Wavelength  l 

■ 0.5  um,  Ax  « 3.15;  S ^ * 6.53955-02 

ext 

■ 6.53955-02  6(ir)  « 3.49286-03  m-l-sr'^ 


Resolution 

Solid  Angle 

. 

BACKSCATTER  i 

(Meters 1 

(Meter) 2 (Steradian) 

(Dimensionless) 

1.0 

1.37-11 

4.19855-14 

2.0 

tf 

7.36763-14 

3.0 

9t 

9.69655-14  ' 

4.0 

It 

1.13437-13  1 

5.0 

It 

1.24412-13  • 

! 

6.0 

It 

1.30991-13  1 

7.0 

It 

1.34087-13  1 

! 

8.0 

If 

! 

1,34455-13 

1 

9.0 

It 

1.32717-13  j 

10.0 

It 

1.29385-13  1 
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Table  D14.  Intensity  Ratios  for  Distribution  2;  Wavelength  ‘ 

- 0.5  un.  Ax  = 3.15;  3^^^.  = 1.99431-01 

B ■ 1.99431-01  B(ir)  = 8.69341-03  m-^-sr-l  i 

sc  j 


1 

Resolution 

Solid  Angle 

BACKSCATTER  i 

(Meters) 

(Meter) 2 (Steradian) 

(Dimensionless)  I 

1.0 

1.37-11 

7.99258-14 

2.0 

M 

1.07274-13 

3.0 

ft 

1.07984-13 

4.0 

It 

i' 

9.66220-14 

t 

5.0 

tt 

8.10518-14  i 

6.0 

It 

6.52710-14  1 

7.0 

tt 

5.11027-14  j 

8.0 

It 

3.91933-14  i 

9.0 

It 

r 

2.95897-14  , 

10.0 

tt 

2.20635-14 

( 


MMifliiiiilMtt 
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Table  DIS.  Intensity  Ratios  for  Distribution  3;  Wavelength 
a 0.5  um,  Ax  = 3.15;  = 3.90475-01  m"^; 

= 3.90475-01  0(Tr)  = 2.49455-02  m'^-sr"^ 


Resolution  Solid  Angle  BACKSCATTER 

(Meters)  (Meter) 2 (Steradian)  (Dimensionless) 


1.0 


2.0 


3.0 


4.0 


5.0 


6.0 


7.0 


8.0 


9.0 


1.37-11  1.56513-13 

" 1.43357-13 

" 9.84801-14 

" 6.01347-14 

" 3.44250-14 

" 1.89187-14 

" 1.01083-14 

" 5.29063-15 

" 2.72582-15 


10.0 


1.38705-15 
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Table  D16.  Intensity  Ratios  for  Distribution  4;  Wavelength 
* 0.5  pm.  Ax  = 3.15;  ®gxt  ~ 2.27353-01 
0^^  = 2.27353-01  m'l;  B(it)  = 1.44447-02  m'^-sr-l 


Resolution 

(Meters) 

Solid  Angle 
(Meter) 2 (Steradian) 

BACKSCATTER 

(Dimensionless) 

1.0 

1.37-11 

1.25589-13 

2.0 

11 

1.59407-13 

3.0 

It 

1.51747-13 

4.0 

II 

1.28405-13 

S.O 

II 

1.01863-13 

6.0 

II 

7.75750-14 

7.0 

II 

5.74371-14 

8.0 

II 

4.16590-14 

9.0 

tf 

2.97430-14 

10.0 

If 

2.09733-14 
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